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I.  Introduction 

This  annual  report  contains  investigations  on  the  properties  of 
nonideal  plasmas,  which  were  carried  through  in  the  period  from  1  November 
1978  to  1  November  1979,  under  ONk  Contract  N0001 4-79-0073 .  Progress 
was  made  in  the  evaluation  of  the  electrical  conductivity  and  the 
thermodynamic  functions  of  nonideal  plasmas.  Research  which  was  pub¬ 
lished  in  this  period  is  not  included. 

Chapter  II  gives  a  dimensional  analysis  of  the  possible  electrical 
conductivities  of  nonideal  classical  and  quantum  plasmas,  with  and 

without  thermal  effects.  Quantum  effects  are  considered  since  the 

21  -3 

electrons  are  noticeably  degenerate  for  densities  n  >  10  cm  and 
temperatures  T  *  10^°K.  In  strongly  nonideal  plasmas,  y  -  7.e? n  ^/KT  >>  1, 
the  thermal  energy  is  negligible  compared  with  the  energy  of  the 
Coulomb  microfields,  i.e.  such  plasmas  behave  like  zero-temperature 
systems . 

Chapter  III  presents  a  kinetic  theory  of  the  electrical  conductivity 

of  nonideal  plasmas,  when  the  electrons  have  i)  a  Maxwell  and  ii)  a 

Fermi  distribution  of  velocities.  This  work  is  based  on  the  quantum 

mechanical  scattering  cross  section  for  an  effective,  shielded 

Coulomb  potential,  which  is  applicable  to  intermediate  nonideal  conditions, 

0.1  <  y  t  10.  For  moderate  degeneracy,  n  <  n,  the  conductivity 

depends  on  a  Coulomb  logarithm,  but  not  for  complete  degeneracy, 

3/  ]  5  i/ 

n>>h,  where  n  =  2 (2nmKT/h2 )  ■-  =  A.  828  x!0  T  ■  . 

Chapter  IV  is  concerned  with  a  statistical  theory  of  Lite  free 
energy  of  nonideal  gaseous  (0  y  <  1)  and  quasi-liquid  (1  <  y  *"  y  ) 
plasmas,  where  yc  is  the  critical  interaction  parameter  of  solid 
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plasma.  The  physical  model  used  considers  short  and  long  range  Coulomb 
interactions  through  quasi-lattice  interactions  and  collective  electron 
and  ion  waves.  The  degeneracy  effects  on  the  free  energy  are  discussed 
for  large  interaction  parametersy. 

Chapter  V  derives  the  Hamilton  function  and  the  canonical  field 
variables  for  a  many-component  plasma  continuum  with  longitudinal 
Coulomb  field  interactions.  By  Fourier  analyzing  the  canonical 
fields  for  random  fluctuations  the  Hamilton  function  can  be  used  to 
calculate  the  free  interaction  energy  of  nonideal  plasmas.  This 
method  has, however ,  not  yet  lead  to  concrete  statistical  applications 
because  of  the  mathematical  difficulties  associated  with  the 
evaluation  of  the  complex  integrals  in  the  multi -dimensional  phase 
space  of  the  Fourier  amplitudes  of  the  canonical  fields. 

Chapter  VI  gives  a  simple  application  of  the  theoretical  approach 
in  Chapter  V  concerned  with  the  evaluation  of  the  distribution  function 
of  the  velocity  fluctuations  in  a  neutral  one-component  gas.  The 
theoretical  distribution  function  is  shown  to  be  in  agreement  with 
the  observations  in  turbulent  gases. 

Chapter  VI 1  is  an  appendix,  in  which  an  unrelated  subject  is 
discussed,  namely  the  propagation  of  stress  relaxation  waves  and  their 
use  for  signal  propagation  and  system  detection  in  water. 

The  main  results  obtained  are  summarized  in  the  abstracts  of 
Chapters  II  -  VII.  Other  research  concerned  with  calculations  of 
the  electrical  conductivity  1)  using  quantum-field  theoretical  methods 
of  solid  state  physics  and  2)  considering  the  effects  of  the  fluctua¬ 
tions  of  the  Coulomb  microfields  on  the  current  transport  could  not  be 
written  up  in  time,  and  will  be  communicated  separetelv.  The  same  holds 
for  tin;  Ph.D.  thesis  of  Mr.  A.H.  Khalfnoui  on  nonideal  plasmas. 
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II.  Dimensional  analysis  of  electrical  conductivity 

OF  NON  IDEAL  CLASSICAL  AND  QUANTUM  PLASMAS 


Abstract 


By  moans  of  dimensional  analysis,  novel  formulas  for  the  electrical 
conductivity  of  nonideal  i)  classical  and  ii)  quantum  plasmas  are  derived 
based  on  the  axioms  of  Dupre"!  In  the  general  case  of  a  nonideal  plasma 
with  partially  degenerate  electrons,  the  conductivity  is  of  the  form  a  * 
Cq  >T  ^  y  B  en  ^  /ra  ^2  ,  where  Cq  is  a  dimensionless  constant,  A  and  B 
are  powers,  and  yT  *  e2n  ^/KT  and  Yq  *  e2n  ^ / (fi2 /ran  ^ )  are  the  reduced 
(7!=1)  nonideality  parameters  of  the  classical  and  quantum  plasma, 
respectively  (e,  m,  n  are  the  charge,  mass,  and  density  of  the  electrons, 
KT  is  the  thermal  energy,  and  it  is  Planck's  constant).  The  known 


conductivities  are  obtained  as  special  cases  of  this  conductivity 

3/  1/ 

rormn in ,  e.g.,  the  conductivity  a  ■  C  (KT)2/m/2e2  of  the  ideal,  classical 

o 

plasma, or  t  lie  conductivity  a  *  C^(me6/fi3)n  ^/KT  of  the  high-temperature 


me t  a 1 . 
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INTRODUCTION 


In  recent  years,  numerous  measurements  of  the  electrical  conductivity 

1—8) 

of  non-ideal  plasmas  have  been  reported, -  which  were  produced  by  oven 

heating  ^ (medium  pressures)  and  shock  wave  heating  — ^  (high  pressures)  of 
alkali  vapors  and  noble  gases,  with  typical  pressures  ranging  from  10* 
to  10~*  bars.  In  spite  of  the  availability  of  an  approximate  kinetic 
equation  for  nonideal  plasmas  which  considers  spatial  and  temporal 
correlations  in  the  collision  operator,  satisfactory  theoretical  explana¬ 
tions  of  the  experimental  conductivity  data  on  nonideal  plasmas  are 
missing  to  date.  The  degree  of  the  nonideality  of  a  fully  ionized  plasma 

is  measured  in  terms  of  the  (dimensionless)  interaction  parameter  y,  which 

,  1/, 

represents  the  ratio  of  average  Coulomb  interaction  (Zezn  J)  and  thermal 
(KT)  energies  (n  =  electron  density,  Z  =  ion  charge  number,  e  =  elementary 
charge) , 

y  =  Ze2n  '  'KT  . 


-3  -1 

In  cgs-unit,  y  =  1.670  x  10  Zn  JT  .  The  conductivity  theory  of  ideal, 

12) 

fully  ionized  plasmas—  agrees  with  the  experimental  data  only  if  y  <<  1. 
For  moderately,  0.1  <  y  <  1,  nonideal  plasmas,  the  ideal  conductivity 
theory  yields  much  too  large  conductivity  values-— -  (KT)  ^2/m  ^2e2Z  fn  A, 
where  l\  =  { l+(D/pQ)2  =D/pq  for  D  >>  pQ ,  p^  =  Ze?/2KT  is  the  impact 
parameter  for  90°  deflections,  and  D  is  the  maximum  impact  parameter. 

The  electric  shielding  length  of  Debye  is  related  to  the  interaction 
parameter  y  by 

D  =  [Z/4n(l+Z)]  ^2y  ^2n 


The  ideal  conductivity  theory  breaks  dc'm  at  higher  electron  densities 
because  tiie  Debye  radius  D  loses  its  physical  meaning  as  an  electric 
shielding  length  and  upper  impact  parameter  when  the  number  of  electrons 


in  the  Debye  sphere,  =  4nD3n/3,  is  no  longer  large  compared  with  one. 
and  y  are  related  by 

SID  =  (4tt/3)  [Z/4tt  (1+Z)  ]  ^  y  ^ . 

For  typical  nonideal  conditions,  n  >  102°cm  3  and  T  =  104°K,  the  Debye 
_8 

radius  is  D  <10  cm,  i.e.,  is  smaller  than  the  atomic  diameter,  which 
shows  that  the  electric  shielding  concept  is  not  applicable  to  proper 
nonideal  plasmas,  y  >  1.  Another  reason  for  the  inapplicability  of  the 
ideal  conductivity  theory  to  proper  nonideal  plasmas  is  its  assumption 
of  successive,  small  binary  interactions,  whereas  in  reality  a  conduction 
electron  experiences  many-body  interactions  for  y  1. 

In  the  following,  we  apply  dimensional  theory  to  the  derivation  of 
new  formulas  for  the  electrical  conductivity  of  (non-relativistic)  ideal 
and  nonideal,  classical  and  quantum  plasmas.  In  the  most  general  case 
of  an  electron  plasma,  the  electrical  conductivity  a  is  a  power  function 
of  the  characteristic  plasma  parameters,  a  ~  ePmCJtir(KT)  Snt  (tf  is  Planck's 
constant).  As  special  cases,  the  conductivity  formulas  for  the  ideal, 
classical,  fully  ionized  plasma  and  the  partially  degenerate,  solid  metal 
are  obtained.  The  derived  formulas  for  nonideal  classical  and  quantum 
plasmas  indicate  the  dependence  of  o  on  y,  which  can  be  compared  with 
the  experimental  observations. 
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THEORETICAL  FOUNDATIONS 

In  a  system  of  reference  in  which  magnetic  fields  are  absent,  a  linear 

4  -4 

electric  current  response  j  =oE  exists,  provLded  that  the  generating 
electric  field  E  is  sufficiently  weak.  For  any  gaseous,  liquid,  or  solid 
medium,  the  electrical  conductivity  o  =  |j|/|F,j  is  given  by 

a  =  (nc2/m)x  (1) 

where  e  is  the  charge,  m  is  the  mass,  n  is  the  density,  and  t  is  the 
(average)  momentum  relaxation  time  of  the  current  carriers.  Because  of 
the  large  ion  mass  m^  >>  m,  the  main  current  carriers  m  a  plasma  are 
the  electrons.  Eq.  (1)  holds  for  any  perturbed  Maxwell  or  Fermi  distribution 
of  the  electrons. 

•  13)  „ 

Dimensional  analysis  is  based  on  the  axioms  of  Dupre.  -  By  axiom 
1),  absolute  numerical  equality  of  quantities  a  and  b  may  exist  only  when 
the  quantities  are  similar  qualitatively.  That  is,  a  general  relation 
may  exist  between  two  quantities  a  and  b  only  when  the  two  quantities 
have  the  same  dimension.  By  axiom  2),  the  ratio  of  the  magnitudes  of 
two  like  quantities  a  and  b  is  independent  of  the  units  used  in  their 
measurement,  provided  that  the  same  units  are  used  for  evaluating  each. 

In  general,  any  measurable  quantity  o(the  secondary  quantity) 
can  be  expressed  in  terms  of  tiiose  appropriate  quantities  a^,  i  =  1,2,...  M 
(the  primary  quantities),  which  affect  the  magnitude  of  o.  The  general 
relationship  between  the  magnitude  of  the  secondary  quantity  a  and  the 
magnitudes  of  the  primary  quantities  a^  is  a  function  of  the  M  arguments 
of  the  form 

o  =  f ( a  j ,  a2,  a3,...aM) 


(2) 
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Application  of  the  axioms  1)  and  2)  to  Eq.  (2)  demonstrates  that  the 


functional  relation  fCa^  is  the  power  function- 
Nj  N2  N3  N  . 

0  =  C  a  *a  za„  . . .  a,  M 
0  1  2  3  :l 


1A) 


(3) 


C  is  a  dimensionless  coefficient,  which  depends  on  the  nature  of  the 
0 

physical  quantity  0,  and  can  only  be  determined  by  means  of  a  detailed 
physical  model.  In  many  cases,  the  order-of-magnitude  of  is  one, 

C  '1. 

o 

In  the  most  general  nonrelativistic  case  of  a  thermal  quantum 
plasma,  of  which  the  classical  thermal  plasma  is  a  special  case,  the 
secondary  conductivity  quantity  0  depends  on  the  dimensional  primary 

quantities  a[  =  e(clectron  charge),  a;>  =  m(electron  mass),  a3  =  ffU’lanck 1  s 
constant ) ,  a[(  =  n(electron  density),  and  a,.  =  KT(t!iermal  energy). 

The  dimensionless  constant  CQ  is  in  general  a  function  of  the  dimensionles 
parameters  p/p^  of  the  plasma, 

C  =  C  fZ.,  m/m. , . . .p/p. , . . . )  .  (A) 

0  01  1  1 

E.g.,Z  is  the  ratio  of  the  magnitudes  of  the  ion  and  electron  charges, 
m/nu  is  the  electron  to  ion  mass  ratio,  etc.  The  electrical  conductivity 
0  and  its  primary  quantities  have  the  following  dimensions  P(L  =  dimension 
of  length,  T  =  dimension  of  time,  M  =  dimension  of  mass)  and  units  U  in 
the  cgs-system: 


Via] 

=  T-1 

♦ 

U[a] 

-1 

=  sec 

V[e] 

= 

* 

U[e] 

3/? 

=  cm  ' gr 

V(m ] 

=  M 

9 

U[m] 

=  gr 

W 

=  ml2t_1 

9 

m] 

=  gr  cm2 

P[n] 

=  IT3 

9 

U[n] 

=  cm 

P[K") 

9  -2 

=  ml2t 

9 

U  [  KT  ] 

-  gr  cm2 

(5) 
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It  is  recognized  that  the  secondary  and  primary  quantities  depend  only 
on  three  basic  dimensions  which  are  independent,  namely  L,  T,  and  M.  ! or 
this  reason,  dimensional  analysis  provides  at  most  three  independent 
equations  for  the  determination  of  the  powers  ,  i  =  1,2,3,.. .M  in  Eq .  (3). 
That  is,  at  most  three  powers  can  be  calculated  while  at  least  M-3  powers 
have  to  be  determined  by  comparison  with  experiments  or  by  physical 
arguments. 

Two  of  the  primary  quantities  are  variables,  namely  n  and  KT 

J  ry  Q 

(K  =•  1.380  x  10  gr  cm" sec  /  K) ,  whereas  the  remaining  three  primary 

-io  Vo  Vo  -l 

quantities  are  elementary  constants  (e  =  4.803  x  10  cm  ^  gr  sec  , 
m  =  9.109  x  10  8gr,  tl  =  1.054  x  10  2?gr  ctn2sec  ). 

The  fundamental  equation  (3)  is  applied  below  to  the  determination 
of  the  electrical  conductivity  of  (nonrelativistic)  non  ideal,  classical 
and  quantum  plasmas,  in  which  the  electrons  are  responsible  for  the 
electric  current  transport.  To  illustrate  the  results,  they  will  be 
expressed  in  terms  of  the  interaction  parameter  y  for  Z  =  1  and  the  order- 


of-magnitude  of  fundamental  electron  energies: 


yt  =  Ec/E?  =  ne2 /KT  , 

E  j  =  KT 
Ec  =  e2nV3 

_2  / 

=  fi2/mn  3 

is  the  order-of-magnitude  of  the  quantum  potential  energy  Q  = 

,  ,  ..  1 L 

-(Tv  /2m)V  (,  z /()  of  the  electron  in  the  plasma.  E  and  E^,  are  the 


(b) 

(7) 

(8) 
(9) 


Coulomb  interaction  and  thermal  energies. 


CLASSICAL  T=0  PLASMA 


In  a  classical.  (fi>-o)  electron  plasma,  in  which  the  thermal  energy  KT 
is  negligible  compared  with  the  Coulomb  interaction  energy,  the  conductivity 
depends  on  the  dimensional  parameters  e,  m,  and  n.  By  Eq .  (3), 


Nj  N2 

a  =  C  e  *m  n 

o 


Hence , 


T'1  -  , 


(10) 


(11) 


(3/2 ) N ,  -  IN,  =  0,  ( 1 / 2 ) N !  +  N2  =  0,  -Ni  =  -l. 

These  are  three  independent  equations,  which  determine  the  powers 
Nj,  N0 ,  and  N3  uniquely, 


(12) 


N,  =  1,  N  =  -V.  , 


Ni  =  *  ■ 


(13) 


(14) 

}h 


By  Eqs.  (10)  and  (i  3)  ,  the  conductivity  of  the  zero-temperature, 

classical  plasma  is 

o  =  C  en  ^  /m  ^ 
a 

Eq.  (14)  indicates  that  a  =  Ca<i^/(4Tr)  ^  increases  proportional  to  n  /2  , 
1/ 

wiiere  w  =  (.4Trne:’/m)  ■’  is  the  plasma  frequency.  This  result  was  first 
P 

derived  by  Euneman—  and  latpr  by  Hamberger  and  Friedman— '  for  an 

electrostatically  turbulent  T  =  0  plasma  by  means  of  semi-quantitative 

1/  1/ 

physical  arguments,  which  give  C^  =  ;r  '  (in^/m)  for  -  I  . 


i 
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CLASSICAL  T>0  PLASMA 

In  a  classical  (h >o)  ,  thermal  (KT>0)  electron  plasma,  the  con- 
ductuvity  depends  on  the  dimensional  parameters  o,  in,  n,  and  KT. 

By  Fq.  (3) 


o  =  C  eNlmN;ynN^  (KT)N|' 

0 


Hence  t 


T~ 1  -  (lVV‘)%N^*3(„|/T-:')S|’  , 


l  .e.  , 


JN]  -  3N3  +  2N4  =  0,  ^-Nj  +  N2  +  N i,  =  0,  -  N,  -  2N,(  =  -1  . 

These  are  three  independent  equations,  which  determine  three  of  the 
powers  in  terms  of  the  fourth. 


Nj  =  1  -  2N,  N2  =  -  }•  ,  N3  = 


_  1  _  1 


N,  N,,  =  N 


(15) 

(16) 

07) 

four 

(18) 


According  to  Eqs.  (15)  and  (18),  the  conductivity  of  the  classical, 
thermal  plasma  is 

«  -  Co  o'-2".,-*  „  5  -  5"  (KT)"  .  (19) 

Collecting  of  powers  of  N  reveals  the  dimensionless  group  contained  in 
Rq .  (19),  and  condenses  the  conductivity  formula  to 

a  =  C(  (KT/e2n  ^3)^  en^/m^2  .  (20) 

For  an  ideal  (classical)  plasma,  o  =  (nc2 /m)t  cannot  depend  on  n 
since  x  ^  n  1  for  binary  e-i  collisions.  Hence,  a  and  N  are  for  an 
ideal  T  •  0  plasma 

=  r  ti rrt  h  h 


a  =  C  (KT)  7  /e2m  ■  ,  N  =  3/2  , 


(21) 


]  2) 

in  agreement  with  kinetic  theory—-  ,  which  shows  that  C 
1/  —  1 

3/4  ( 2 xr )  ^ ZCnA  ~  10  .  For  N  =  0,  Rq .  (20)  reduces  to  o  of  the  classical 


T  =  0  plasma,  Rq.  (14). 


For  a  nonideal  (classical)  plasma  with  v-body  interactions,  we 
have  t  “  n1  V  and  a  -  (ne^/m)r  *  n?  V,  with  v  >  2,  i.e.  N  =  3v  -  9/2  -•  1/2. 


Thus,  we  find  for  the  conductivity  of  nonideal,  classical  plasmas 

„  -N  V?  ,  V?  ... 

o  =  en  //m  '  ,  N  ■  3/2  , 

where  Y,r=  E../E  t'le  nonideality  parameter  defined  in  Eqs.  (6)  - 
l  1 


(22) 


(S) .  F.q.  (22)  expresses  the  important  result  that  the  conductivity  of 

-N  ty> 

a  nonideal,  classical  plasma,  o  =  / 2rt  2  ,  decreases  proportional  to 

-N 

Y^-  witii  increasing  nonideality  since  N  >  3/2.  The  exact  value  of  N 
can  bi  determined  by  comparison  with  experimental  data. 


1, 


T  =  0  QUANTUM  plasma 

In  a  completely  degenerate  electron  plasma,  F^.  •<  Eq,  the  conductivity 
depends  on  the  dimensional  parameters  e,  m,  n,  and  fi,  but  not  on  KT.  By 
Eq.  (3), 


„  N,  No  N,.N4 
o  =  C  e  1  m  ^  n  •  ft  H 


Hence , 


l  .  e .  , 

jN [  -  3N3  +  2N4  =  0,  jNj  +  N2  +  N4  =  0,  -N,  -  N4  =  -1 . 

These  are  three  independent  equations,  which  determine  three  of  the 
four  powers  in  terms  of  the  fourth, 

N,  =  1-2N,  N2  =  --  -  N,  N3  =  i  +  In,  N4  -  2N. 

Substitution  of  Eq.  (26)  into  (23)  yields  as  conductivity  of  the 
completely  degenerate  electron  plasma 


a  =  C 


v* 


For  N  =  3/2,  Eq.  (27)  leads  to  the  conductivity  of  a  solid  metal 


at  T  =  0, 


=  C  tf^n/e?m^  , 


N  =  3/2  , 


17)  -I 

where  •  C  <*  Z  .  For  N=0,  Eq.  (27)  reduces  to  a  of  the  classical 


T=0  plasma,  Eq.  (14). 

For  complete  degeneracy,  E^,  is  negligible  compared  with  E^  ~  -mVp  , 
1/ 

where  Up  -  (tf/m)n  *  is  the  Fermi  velocity.  For  this  reason,  Eq .  (27) 
is  rewritten  in  the  form 

o  =  C(j  en  ^/m  ^  (29 


where 


Yq  =  K(./Kq  =  e?n  ^ / (1v  / mn  ^ ) 


(30) 
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is  the  nonideality  parameter  of  the  completely  degenerate  electron  plasma 
for  Z  =  1.  Since  N  =  3/2  for  the  T=0  metal,  it  is  to  be  expected  that 
N>0  for  the  completely  degenerate  electron  plasma,  i.e.  its  conductivity 
decreases  with  increasing  nonideality  y  ,  Eq.  (29).  This  formula  is 
useful  for  the  interpretation  of  conductivity  data  of  completely 
degenerate  electron  plasmas,  with  N  as  adjustable  parameter. 
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T  -0  QUANTUM  PLASMA 


In  a  partially  degenerate  electron  plasma,  E,j,  £  E^ ,  the  conductivity 

depends  on  the  dimensional  parameters  e,  m,  n,  KT,  and  fi.  By  Eq.  (3), 

„  N,  N  '  N3  N4i.N5 

o^C^e  '-n  J(Kl)  4fi  J  .  (31) 

Hence , 

T~  *  =  (L^M1/2T-1)NlMN2L~,N3(ML2T_?)N'' (ML?T_))N,i  ,  (32) 

i  .  e  .  , 


-N! 

-  3N3  + 

2N4 

+ 

2Ni;  =  0 

9 

+  N2  + 

n4 

+ 

N5  =  0 

9 

(33) 

-N, 

-  2N„  - 

n5 

= 

-1 

. 

These  are  three  independent  equations,  which  permit  to  express  three 
of  the  five  powers  in  terms  of  the  remaining  powers, 

Nj  =  1  -  2A  -  2B  ,  N2  =  -  i  -  B,  N3  =  5  -  |a  +  |b, 

N4  A,  N5  =  2B  .  (34) 

Combining  of  Eq.  (31)  with  (34)  gives  for  the  conductivity  of  the 
T>0  quantum  plasma 


■=•  w  (^)B 


(33) 


whore  A  and  B  are  powers  which  can  not  be  determined  by  dimensional 
reasoning . 

For  A  =  -1  and  B  =  -3/2,  Eq.  (35)  yields  the  conductivity  of  solid 


metals  at  temperatures  T>0, 
me6  n1/s 

a  =  C  — - —  ,  A  =  -1,  B  =  -3/2, 


0  *3 


(36) 


KT 


where— a  Z  ^ .  Eq.  (36)  expresses  the  1/T  -  lav^^of  the 
metallic  conductivity  at  "high  temperatures". 

Eq.  (35)  contains  the  dimensionless  groups  and  Yq»  which  permit 
to  rewrite  the  conductivity  formula  as 

A  -B  V2/  V2 

en  /m  . 


0  =  Co  YT  YQ 


(37) 


Since  A  =  -1  for  T>0  metals  [Eq.  (36)]  and  B  =  3/2  for  T=0  metals[Eq.  (28)] 
one  can  speculate  that  A  >  -1  and  B>0  for  nonideal,  partially  degenerate* 
plasmas.  A  and  B  can  readily  be  determined  by  means  of  conductivity 
measurements  for  nonideal  quantum  plasmas.  The  theoretical  determination 
of  the  nonlinear  dependence  of  o  on  y^  and  Yq  in  Eq.  (37)  from  a 
physical  model  of  more  than  two-body  interactions  is  left  to  future 
research . 

In  conclusion,  it  is  noted  that  we  have  derived  new  conductivity 
formulas  for  nonideal  classical  T>0  plasmas  [Eq.  (20)  or  (22)],  completely 
degenerate  plasmas  [Eq.  (27)  or  (29)],  and  nonideal  T>0  quantum  plasmas 
[Eq.  (35)  or  (37)].  These  formulas  can  be  used  to  interprete  conductivity 
measurements  on  nonideal  plasmas.  Once  the  still  undetermined  powers 
N,  A,  and  B  are  known  empirically,  it  should  be  possible  to  develop  a 
conductivity  theory  for  nonideal  plasmas  which  provides  an  explanation 
of  the  y  and  y.  dependence  from  first  principles. 
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III,  Electrical  conductivity  of  nonideal  plasma 


ABSTRACT 

The  electrical  conductivity  of  fully  ionized,  moderately  nonideal.  plasmas 

2  173 

with  interaction  parameters  0.1  <  y  S  1,  where  y  ■  Ze  n  /KT  is  the  ratio 

of  Coulomb  and  thermal  energies,  is  calculated  for  displaced  Maxwell  and 

Fermi  electron  distributions,  respectively.  The  electrons  are  scattered  by 

-1  -1/3 

an  effective  Coulomb  potential  4>(r)  *  Zer  exp(-r/<$),  with  6  *  (3n/4nZ)  ' 

the  mean  ion  distance,  which  considers  binary  (0 <  r  $  6)  and  many-body 
(<$  <  r  <  »)  interactions.  It  is  shown  that  the  resulting  conductivity  formula 
is  applicable  to  densities  up  to  four  orders  of  magnitude  higher  than  those 
of  the  ideal  conductivity  theory,  which  breaks  down  at  higher  densities 
because  the  Debye  radius  loses  its  physical  meaning  as  a  shielding  length  and 
upper  impact  parameter.  In  the  limit  of  complete  degeneracy,  the  conductivity 
formula  reduces  to  that  of  a  solid  metal. 


18 


.1-3) 


lations^  give  an  electrical  conductivity  of  the  form  a  ~  (KT)  ^/^/m^/^e‘'Z('ttA 


INTRODUCTION 

The  theory  of  the  electrical  conductivity  of  fully  ionized  plasmas'  J/  based 
on  the  Boltzmann  equation,  the  Fokker-Planek  equation  (derived  by  expanding  the 
binary  collision  integral  for  Lhe  small,  successive  velocity  changes  of 
Coulomb  scattering),  or  the  Lenard-Ba 1 escu  equation  (taking  into  account  the 
dielectric  properties  of  the  medium)  is  in  agreement  with  the  experimental  data 
for  rarefied  high-temperaLure  plasmas,  y<<1.  The  interaction  parameter  is 
defined  as  the  ratio  of  (average)  Coulomb  interaction  (Ze^n^^)  and  thermal  (id) 

energies  (n  is  the  electron  density  and  Z  the  ion  charge  number), 

2  1/3  -3  1/3 

y  =  Zc n '  /KT  =  1.670  x  10  Zn  /T 

in  cgs-units  which  will  be  used  throughout.  The  conventional  transport  calcu- 

i  form  a 

2  1/2 

for  classical  ideal  plasmas,  where  A  =  [1  +  (D/pn)  ]  ~  l)/p  for  D>^po.  D  is 

2  1/2 

the  maximum  impact  parameter  (Debye  length),  D  =  fKT/4ire  (1  +  Z)n]  ,  and  p 

2 

Is  the  average  impact  parameter  for  90°  deflections  (Landau  length),  p  =  Ze  /2KT. 

The  condition,  A>>1  or  £nA~1.0^  is  satisfied  only  for  not  too  low  temperatures 

4) 

T  and  not  too  high  densities  n.  Conductivity  formulas  with  this  Coulomb 
logarithm  break/down  for  large  interaction  parameters  y  and  densities  n,  since 
Lhe  Debye  radius 

l>-  |Z/4.(1  W.)l1/V'/Vl/3 

“8 

becomes  smaller  than  the  atomic  dimension  10  cm  and,  thus,  completely  loses 
its  physical  meaning  as  an  electric  shielding  length  and  maximum  impact  parameter. 
E.g.,  tor  T  =  10^  °K,  y>-10^and  D<10  ^cm  if  n>10^cm  Moderately  nonideal 
plasmas  with  y~l  are  readily  generated  through  shock  wave  compression  and  exhibit 
conductivities  of  the  order  u  10^  -  10^mho/cm^  which  are  much  smaller  than 
those  which  would  be  obtained  by  applying  the  conductivity  formula  for  ideal 
plasmas  in  the  nonideal  regime. 
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Although  there  are  some  bulk  measurements  of  the  eleetric.il  conductivity 

5—8} 

of  nonideal  cesium  and  noble  gas  plasmas  available  ,  theoretical  explanations 

of  these  results  are  still  missing.  The  momentum  and  ncrgy  transport  in 

9) 

weakly  nonideal  plasmas, y  «  1,  was  treated  by  Wilhelm  by  means  of  an 
experimentally  shielded  Coulomb  potential,  which  permits  Lo  consider  not  only 
short-range  binary  (r  <  D)  but  also  long-range  many-hody  (r  D)  interactions.  This 
interaction  model  ;a s  used  shortly  afterwards  by  Rogov  for  the  calculation 
of  the  conductivity  of  weakly,  nonideal  argon  and  xenon  plasmas  with  Debye 
shielding. 

For  moderately  nonideal  plasmas,  0.1  <y  S  1,  various  phenomenological  approaches 

have  been  used  to  extend  the  conductivity  formula  of  ideal  plasmas,  e.g.,  Goldbach 

et  al  multiply  the  Debye  length  D  with  a  free  parameter  x(p)which  is  chosen 

to  match  the  experimental  data,  i.e.  to  compensate  for  the  too  rapid  decrease  of 

D  with  pressure.  A  kinetic  equation  has  been  proposed  for  nonideal  plasmas 

12) 

by  Klimontovich  ,  which  considers  spatial  correlations  and  temporal  retardation 
in  the  collision  integrals.  This  equation  has  not  yet  lead  to  transport 
coefficients  because  of  the  mathematical  difficulties  associated  with  its  solution. 

In  the  following,  the  momentum  relaxation  time  and  the  electrical  conductivity 
of  (i)  classical  and  (ii)  degenerate  plasmas  is  calculated  for  intermediate  non¬ 
ideal  condi tions,  0.1  <  >  S  1.  For  this  region  of  interaction,  the  concept  of 


Debye  shielding  already  breaks  down  since  the  number  of  particles  in  the  Debye 

3 

sphere  4irD  /3  is  no  longer  large  compared  with  one  for  y  >  0.1.  This  difficulty 
can  not  be  remidied  by  replacing  D  with  the  quantum  mechanical  shielding  length 


(even  when  thermal  effects  are  included)’ 


Dp  =  (nao/kF) 


a  -Sly,  k, 
o  mo2  1'  8,1 


which  is  of  the  same  order  as  D  in  most  high  pressure  plasmas,  e.g.  D  ~  10  cm 

r 

for  n  =  lO^  cm  From  the  definition  of  the  mean  particle  distance,  it  is  clear 


that  the  mean  ion  distance  6  ~  n  separates  the  region  in  which  an  electron 
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experiences  few  -  body  encounters  (r  <  ( )  from  the  region  in  which  an  electron 

experiences  many  -  body  interactions  (r  >  5)  in  a  nonidea]  plasmas,  as  long  as 
—8  24  —  3 

5  >  10  cm  (n^  <  10“  cm  ).  Thus,  the  mean  ion  distance  evolves  naturally 
as  the  characteristic  interaction  distance  for  nonideal  plasmas,  for  which 
Debye  and  Fermi  shielding  fail. 

We  calculate  first  the  electrical  conductivity  of  plasmas  with  (i)  Maxwell 

and  (ii)  Fermi  distributions  of  the  electrons,  when  all  ions  have  the  same 

charge  number.  Then,  we  generalize  the  conductivity  formulas  for  plasmas 

with  several  ion  components  of  charge  number  Z^,  i  =  1,  2,  ...  .  The  electrons 

are  assumed  to  be  scattered  by  the  exponentially  decaying  Coulomb  potential 

p  =  Zer  ^  exp(-r/<5)  which  takes  many  -  body  interactions  at  distances  r  >  6 

into  account.  The  considerations  are  applicable  only  to  moderately  nonideal 

24 

conditions,  0.1  <  y  <  1,  up  to  densities  n  <<  10  cm. 

24  -3 

For  solid  state  densities  and  larger  densities,  n  >  10  cm  ,  we  have 
-8 

6  <  10  cm  and  the  chosen  interaction  potential  t(r)  is  no  longer  valid.  For 

24  -3 

densities  n  >10  cm  ,  the  Coulomb  field  of  the  ions  is  changed  by  polari¬ 
zation  of  the  bound  electron  cloud,  so  that  the  free  electrons  are  scattered 
as  in  solids  by  the  ions.  For  this  reason,  plasmas  in  the  solid  phase  are 
not  treated  herein. 
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PHYSICAL  FOUNDATIONS 

The  electrical  conductivity  o  of  ,-my  gaseous,  liquid,  or  solid  medium,  in 

which  the  electrical  current  transport  is  due  to  electrons,  is  proportional 

to  the  electron  density  n  and  the  relaxation  time  t  of  the  average  momentum 

<mv  of  the  electrons  (m  is  the  electron  mass  and  e  >  0  is  the  elementary 
o  J 

charge) 

2 

o  =  (ne  /nt)  x  (1) 

hq.(l)  reduces  the  calculation  of  a  to  the  evaluation  of  i .  The  relaxation  time 
r  is  determined  by  Lhe  scattering  potential  and  the  (classical  or  quantum 
statistical)  kinetics  of  the  electron  gas  in  the  electric  field. 

In  proper  nonideal  plasmas,  the  region  0  <  r  <  6  of  binary  and  few  -  body 
collisons  and  the  region  6  <  r  <  «  of  many  -  body  interactions  are  bounded  by 
the  mean  ion  radius, 

(4n53/3)n.  =1,  6  =  (3Z/4tm)l/3,  (2) 

since  electric  Debye  shielding  exists  only  for  weakly  nonideal  conditions,  y  <<  1 . 
For  this  reason,  the  effective  Coulomb  potential  of  Z  times  charged  ions  is  in 
plasmas  of  intermediate  nonideality 

<1> ( r )  =  Zer  ^  exp(-r/6),  0  <  r  <  <®  ,  0.1  <  y  <  1.  (3) 

Eq.(3)  is  no  longer  applicable  to  plasmas  with  densities  comparable  to  solids.  It 

contains  the  binary  and  few  -  body  collisions  at  distances  0  <  r  <  5  and  the 

many  -  body  interactions  at  distances  6  <  r  < 

The  differential  cross  section  o(0,g)  for  the  scattering  (g  -*•  g  *)  of 

14) 

electrons  by  the  potential  (3)  is  in  the  center  of  mass  system 

o(0, g)  =  (Ze2/2m)2/[g2  sin2(0/2)  +  u2]2  (4) 

where 

0  =  Mg,g*),  g  =  vc  -  Vj_,  8*  =  v*  -  v£*  (5) 

and 


u  =  Ti/2m4 


(6) 


Tile  electron  and  ion  velocities  oefore  and  alter  l lie  interaction  are  designated 
-y  -> 

by  v  .  and  v*  .  ,  respectively.  The  speed  u  corresponds  to  a  de  Broglie 

G  )  1  G  t  l 

wave  length  of  the  order  \  -  6.  For  u  >•  0  or  ■  -  Kg. (4)  reduces  to  the 

Rutherford  cross  section. 

The  scattering  cross  section  u(0,g)  is  strictly  valid  only  in  the  Born 
approximation^^.  Contrary  to  what  one  might  expect  in  general  lor  the  latter, 
Eq.(4)  describes  in  good  approximation  the  scattering  in  the  exponential iy 
decaying  potential  (3)  because  of  the  peculiarity  of  the  Coulomb  interaction. 

The  Coulomb  interalion  (J>  -  1/r  has  the  unique  property  that  the  Born  approximation 
and  tlie  exact  wave  mechanical  approach  give  the  same  scattering  cross  section^ 
(identical  with  the  Rutherford  formula).  In  the  region  0  <  r  <  6 ,  the 
interaction  potential  (1)  is  practically  Couiombic,  and  thus  the  Born  approxi¬ 
mation  gives  the  correct  solution.  In  the  region  6  <  r  <  »,  the  interaction 
potential  (1)  is  effectively  screened,  i.e.,  the  Born  approximation  gives  the 
correct  solution  because  ^(r)  is  small.  1 n  the  transition  zone  r  ‘  n,  the 
Born  approximation  bolds  fairly  well  for  reasons  of  continuity. 

The  relaxation  time  i  is  obtained  by  evaluation  of  the  collision  integrals 

-► 

for  tlie  electron  momentum  mv  for  the  (i)  classical  and  ( i  i )  degenerate  plasma, 

respectively.  Both  in  the  cases  of  classical  and  Fermi  statistics,  the  particle 
-y  >■ 

velocities  v  .  and  v*  .  before  and  after  the  interaction  are  interrelated  by  the 
tj *  i-  g  ,  i 

conservation  equations  for  momentum  and  energy. 


CONDUCTIVITY  OF  CLASSICAL  i’LASMA 


According  to  kinetic  theory,  tin-  average  momentum  density  n  m(-'v  ■ 

e  1 

exchanged  per  unit  Lime  between  electrons  and  ions,  interacting,  with  the  Coulomb 
potential  tl),  is  given  by  the  collision  integral  lor  mv  ,  which  determines 

U 

the  momentum  relaxation  time  t  , 

->  •  > 

“Illll(<v  >  -  <V.>)/l  “ 

o  1  ' 

m  /.../  V  |f  (V*)  1  .  (v*)  -  f  (V  )  f .  ( v . )  ] ,;  .(',g)  d:,;  dv  dv.  .  (7) 

J  1  e  e  e  1  l  e  e  i  i  e  i 

The  scattering  cross  section  c  (>,  g)  is  given  in  L-  •  j .  C  4 )  .and  the  sol  iti  angle  element 

is  d'-?  =  sin  0  di)  do.  In  response  to  an  applied  electric  fit  Id  E,  Lite  electrons 
t  ►  >• 

and  ions  dr-'ft  with  velocities  <v  and  <v .  >  so  that  their  distribution  1  suctions 

e  i 


are  displaced  Maxwell ians, 

f  ,('■  ,)  =  n  (m  ,/2-tKT  )  ^~  oxp[-~  m  (v  -  <v  .  •)‘7KT  ) 

S  S  S  S  S  iL  S  !■>  ,  t  s 


■  ,  i . 


(«) 


E<( .  (8)  represents  a  5-moment-approximnt  ion  to  the  nonecpii  1  ibrium  solution  of  Lite 
Boltzmann  equation.'1*'  The  perturbations  of  f  (v  )  due  to  viscous  stresses  and 

S  ii 

heat  flows  are  neglected  in  Eq.(8),  since  they  yield  only  corrections  of 
higher  order  to  the  conductivity. 

9) 

The  collision  integral  (7)  is  integrated  by  standard  methods  for  subsonic 
drift  velocities,  j  <v  >  -  <v  >  |  <  (2KT/ra) *  with  the  usual  approx i mat  ions 


(m  .  =  m  m .  /  (m  +  m  . ) 
e  L  e  l  e  l 


m  :  in,  T  =  m  I  (T  /m  )  +  (T./m.)]  T 
e  e  S  t  s'  e  e  J  .  .. 


T;  . 


For  supersonic  drift  velocities,  a  linear  response  j  =  oL  between  current  d.  nsity 

9 ) 

j  and  electric  field  E  does  no  longer  exist.  The  resulting  relaxation  time  is 


m 


g iven  by: 

-1  8 .. 


"(2KT/iim)  1/2  u.  (> 


^(Ze2/  KT)2  L 


L,  (A  ') 


(h) 

(10) 

OD 


wlieri 


I’KT/mu2  -  <Hiii/'ir  )  (  m/  1Z)  V,  | 


(12) 


is  t  ho  exponential  integral  of  order  one  (T  =  0.477  ...  -  Killer's  constant).  ' 


The  latter  satisfies  the  inequalities,  for  x  0, 


+  2/x)  <  eXE  L  (x)  -  t'»U  +  1/x),  (1  +  x)  '  ,.'i  j(x)  <  X  ( 

According  !y ,  E<| .  ( 1  1 )  gives  formally  lor  sm.il]  and  1  urge  :>  reunion  tu  =  A  K 

L  ■-  tn  A  ,  A  >  1  ,  ( 

r.  -  a  ,  a  «  i  ( 


Combining  of  Eqs.  (9)  -  (12)  with  ]•;  | .  (1)  yields  tlx-  desind  electric 
conductivity  of  the  classical  plasma  of  iniormed  i.ile  nonideality,  (J.  I  <  c  <  1  : 

o  =  3(K.T)-J/2/2(2tim)1/2  e2  El.  ( 

where 

-  •’  _ >  /  a 

h  =  f'n[8mti  ~  (4ttn/3Z)  KT  | ,  A  >•  I,  ( 

-2  -2/3 

1.  =  Binlf  (4irn/3Z)  '  KT  ,  \  «  \.  ( 

by  Eqs .  (15)  -  (16). It  is  instructive  te  rewrite  A  in  terms  of  the  thermal  and 
quantum  potential  energies  of  an  electron, 

A  -  8ET/EQ,  Et  =  KT,  E  =  ir/mh2  .  ( 


The  conductivity  formula  (17) 
plasma^  mainly  through  the  term 
logarithm,  L  =  Cn  A  ,  for  A  ■  •  1,  i 
which  the  plasma  is  nondegenerate, 

A  =  3.482  x  K)1 1 (n/Z)_J/3  T 
The  corresponding  argument  2KT 


differs  from  the  conductivity  of  the  ideal 
L.  The  latter  has  the  form  of  a  Coulomb 
.e.  for  all  densities  n  and  temperatures  T 
E.j,  E  ,  Kq.  (20)  .  Numerically, 

3  1-3) 

D/Z.'-  of  L  lie  ideal  Ciuilomh  logarithm  , 


for 


( 


Table  1  compares  A  of  the  nonideal  plasma  and  of  the  ideal  plasma  for  large 

densities  n  and  the  typical  temperature  i  ;  10  ,JK.  It  is  seen  that  hi  A 

i)  of 

the  ideal  plasma  is  unacceptably  small  for  densities  n  >  I  (J  '8  cm  ^ ,  whereas 

,  22-1 
t  II  '  ol  the  nonideal  plasma  has  reasonable  values  up  to  densities  n  ■:  10  cm  , 

if  T  =  10*  °K.  The  conductivity  formula  (17)  holds,  therefore,  for  densities  n 

up  to  A  orders  of  magnitude  higher  than  the  conduct  ivitv  formula  of  the  ideal 

plasma.  l’.q.(17)  is  not  applicable  to  n  -  T  r.  j-.utvi  for  which  A  [,  L.e. 

E  *<  F.  ,  Eq.  (20)  ,  since  in  this  case  the  i  I  •  t  i  mi;;  would  he  degenerate. 


nlcm  J] 


TABLE  1:  A  and  A  versus  n  f.  r  1'  -  10  °K  and  7.  =  1. 


3.482  x  103  1.616  x  102  0.750  x  10’  0.348  x  10° 

1.035  x  L01  1.035  x  10°  1.035  x  10_1  1.035  x  10~2 


The  conductivity  formula  (17)  becomes  in  cgs  -  units  or  practical  units 


(y  x  10^  see  3  =  1  mho  citt'b , 


•  =  l . '394  x  108  T3^2  /Z  hi  A  [sec  ‘  ]  =  1.549  x  10  4  T 3,/2/Z  hi  A  [mho  cm  ]]  (23) 
where  A  Is  given  in  Eq.(21).  Accordingly,  if  1  =  !0+  °K  and  Z  =  1,  o  =  1.899  x  10 
mho  enf1  for  n  =  10J8  cm”3  and  o  -  2.097  x  101  mho  cm-1  for  n  =  102°  cm  3. 


CONDUCTIVITY  01'  QUANTUM  IM.ASMA 


The  electrons  in  a  plasm.-,  become  degenei  ‘i  .v  if  their  thermal  DeBrogl.  ie 


vave  lenj’.th  is  larger  than  the  mean  electron  I  i  si  .nice,  i.e.  al  densities  (h=2iiti) 


n  •  2  ( Turn  KT 


/h2)3/2  =  4ji28  x  1015  T3/"\ 


2l  ..„-3 


l l.g.,  for  T  =  1 O*4  K,  degeneracy  requires  n  ;  x  \(r  tin  .  In  view  of  their 


arge  mass  m . 


m,  the  ions  can  be  treated  as  classical.  The  momentum  relaxation 


time  i  of  the  degenerate  electron  gas  is  determined  by  the  quantum  statistical 

-►  19) 

collision  integral  for  niv  , 

e 

->  •  -V 

-  nm(<v  -•  -  <v.a)/i  = 


v  {f  (v*)f.(v*)[l  -  f  (v  )]  -  f  (v  )  f .  (v . )  [  I  -  Vf  'v*)l 
J  J  e  e  err  2  e  e  e  e  t  l  2  e  e 


x  g  o (0 , g)  dW  dv^  dv^  ( 

where  the  scattering  cross  section  o^.g)  between  electrons  and  ions  is  given  by 
Eq.  (4) .  The  solutions  to  the  velocity  distributions  are  the  displaced 
Maxwellian  (8)  for  the  ions  (s  =  i)  and  tiie  5-moment  Fermi  approximation  for  the 
electrons , 

f  (v  )  =  2 (m/h)  !  tl  +  exp [Jan(v  -  v  •)"  -  u|/KT  i  1  ( 

e  e  2  e  e 

The  chemical  potential  p  =  p(n,  T)  is  determined  by  the  integral  functional 


n  =  Jf„(c  ;  p)  dc  . 
e  e  e 


Again,  a  linear  response  j  =  oE  exists  for  small  drift  velocities  <v 

e ,  i 

or  weak  electric  fields  E.  integration  of  Eg. (25)  yields,  after  standard 

approximations,  for  the  relaxation  time  of  the  degenerate  electron  gas: 

-1  _  8.2  1/2  e4Z2nj  LQ _ _ 

3  m  (RT)3/2R(n,T) 


where 


.  _  16  m  4  an. -2/3  n 

\)  ~  3  ^2  (~3Z)  Kf  a(n’° 


2  7 


and 


Q(n,T)  =  |(1  +  2  2  "  +  ...),  n  <•  n  ( 1')  ,  (30) 

2-  n 


R(  n,T)  = 

ii 

(1  +  2  2  ~  +  ...).  n  n (T) , 

n 

(31) 

but 

Q(n,T) 

3  3n 1 /2  n. 
*  '  n 

m\i  +  -5r(-34-  ?)■*'’•*  -  .1. 

^  n 

it  '  •  n(T)  , 

(32) 

R(n,T) 

1  3n1/2  n 

0  *  /.  . :  ' 

,  q  2  .  1/2  -4/3 

[  1  +  ~  (--;--  -:) 

n  -  •  n (T)  . 

(33) 

Equations  (30)  -  (31)  and  Eqs.  (32)  -  (33)  result  from  expansions  of  the  Eermi 
distribution  (26)  in  the  collision  integral  (25)  for:  densities  n  «  n(T)  and 
n  ->'■  n(T)  ,  where 


2  3/2 

n  (T)  =  2  (2-nm  KT/h")J/ 


( 34) 


is  tire  critical  density  which  separates  the  degenerate  and  non-degenerate  regimes. 
These  series  are  based  on  expansions  of  the  normalization  integral  of  Eq. (26) , 
which  gives  the  chemical  potential  m  explicitly  as  a  function  of  n  and  T, 

-£=  =  t»l{"  [1  +  2~3/2(")1  +  4  "  3 "J/2)(")2  +...],  n  '  n(T),  (35) 

Kin  n  ^  n 


P_  _  ,  3  ii 
KT  '  A 


0  ...  2  ,  1/2  .  A  1/2  H/„ 

n.  2/3  .  _n  ,3n  _  n.--A/3 _ *  .Jh__  n  -8/3 

1  12  1  4  „  80  A  '  „ 


n  -  n(T)  (36) 


n  ^  n 

Combining  of  the  conductivity  formula  in  Eq.(l)  with  the  relaxation  time  of 

Eq(27)  yields  for  the  electrical  conductivity  of  the  degenerate  electron  plasma 
of  intermediate  nonideality,  0.1  <  Y  :>  1: 


3 (KT) 3/2  R(n,T) 
8(2m)1/2e2Z  L 

where  I.  is  given  by  Eq.(29), 


(37) 


In  the  limiting  cases  of  large  and  small  values  of 


I.,  =  0  (n  ,T)  |  , 

»  >  i 

h()  =  -2A0(n,T), 


A  0 
0 


(38) 


(39) 
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since  A  =  (2/3) AQ  by  comparison  of  Eqs. (29)  and  (12). 

For  n /n  ->  0,  Eqs. (37)  and  (38)  reduce  to  the.  classical  conductivity,  Eqs. 
(1/)  and  (18),  since  R(n,T)  >  4//n  and  0(n,T)  <•  3/2  for  n/n  0  by  Eqs.  (30) 

-  (31).  On  t lie  other  hand,  Eqs.  (37)  and  (3'))  p.ive  in  the  limit  of  coniplete 
degeneracy,  n/n  ■>■ 


0  "9  ~2  2 

2  inn  e  Z  I, 


where 


,  =3  3  L/3Z_ 

Q  40',2ir  3 


lq  ;<1 


by  Eqs.  (32)  and  (33).  Euqations  (40)  -  (41)  combine  to  the  conductivity  lormula 


.15tt2  .2  a  1/3  h^n 

.,6  C3  ’  2  2  _5/3 

2  m  e  Z 


Equation  (42)  agrees  with  the  expression  for  the  conductivity  of  a  low  temperature 

,  20) 

metal. — 

For  numerical  evaluations,  Eq.(42)  is  stated  in  cgs  units  (sec  )  and 


practical  units  (mho  cm  ), 

_  /.  /.no  ia’^  n  -5/3 


LB  „  —3 / 3 


j  -  4.498  x  10  1 '  n  [sec  )  =  4.998  x  10  Z  a  [mho  cm  |. 

Accordingly,  a  4  x  10^  mho  cm-1  for  n  =  10^  cm' J  and  u  5  x  |o  mho  cm  ' 
for  n  =  lO*"4  cm  3?  if  z  =  1. 


2<> 


GENERALIZATION 

Nonideal  plasmas  exhibit  frequently  not  onl\  a  high  degree  of 

single  ionization  but  also  multiple  ionization,  due  to  lowering  of 

the  ionization  energies  by  the  internal  Coulomb  fields,  and  over 

lapping  of  the  atomic  wave  functions  at  sufficiently  high  pressures. 

In  an  electrically  neutral  plasma  with  N  species  of  ions  (i)  of 

charge  Z.e  and  density  n^,  the  electron  density  n  and  entire  ion 

density  n(i)  are  related  by 

N  N 

n  =  l  ^|ni >  n(i)  =  l  n.  .  (44) 

i=l  i=l 

The  characteristic  interaction  radius  6  of  the  Coulomb  field  of  each 
ion  is  within  the  many-component  plasma 

6  =  f  4im(i)/  .  (45) 

Since  the  probabilities  for  interaction  between  the  electrons 


and  ions  of  type  i  =  i,  2,...N  are  additive,  the  momentum  relaxation 
time  of  the  electrons  is 


where  t.  is  the  relaxation  time  of  the  momentum  exchange  of  the 

l  ' 

electrons  with  ions  of  type  i. 


Classical  Plasma.  By  Eqs.  (1),  (46),  and  (9) -(10),  the  conductivity 


of  the  many-component  plasma  with  Maxwell  electron  distribution  is 
in  the  Lorentz  approximation 


2  (2irm)  •  e 


3  (KT) 
-1/—  N- 

I  '«  V 


(nj/n)Zj  •!. 
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I,  is  given  by  Eq.  (11)  in  general,  and  the  approximations  (15)  and 
(16)  in  the  limits  A  >>  1  and  A  <<  1,  respectively,  where  now 

A  =  8KT/ (Ti2  /m62 )  .  (48) 

Quantum  Plasma.  By  Eqs.  (1),  (46),  and  (27),  the  conductivity 
of  the  many- component  plasma  with  Fermi  electron  distribution  is  in 


the  Lorentz  approximation 


3/ 

3 ( KT )  2  R(n ,T) 


!;  N 

8 (2m)  '•  e2  T  (n./n)Z.*L 

i=i  1  1  Q 


Lq  is  given  by  Eq.  (18)  in  general,  and  the  approximations  (38) 


and  (39)  in  the  limits  A^  >-■  1  and  A^ 


.1  ,  respectively,  where 


Aq  =  (16/3)KTO(n,T)/(-ri2/m6’)  .  (50) 

For  a  plasma  with  only  one  ion  component,  we  have  N  =  1, 

7,^  =  Z,  n./n  =  Z  1  ,  and  Eqs.  (47)  and  (49)  reduce  to  the  previous 
formulas  (17)  and  (37),  respectively.  In  case  of  the  many-compnent 


plasma,  not  only  the  electron  density  n  but  also  all  ion  densities  n^ 
have  to  be  known  for  the  evaluation  of  the  conductivity. 
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IV.  Statistical  thermodynamics  of  nonideal  plasma 


Abstract 

A  quantum  statistical  theory  of  the  free  energy  of  a  nonideal  electron-ion 
plasma  is  developed  for  arbitrary  interaction  parameters  0  <  y  <  y 

c 

1/. 

(y  =  Ze' n  J/KT  is  the  ratioof  mean  Coulomb  interaction  and  thermal 

energies),  which  takes  into  account  the  energy  eigenvalues  of  (i) 

the  thermal  translational  particle  motions,  (ii)  the  random  collective 

electron  and  ion  motions,  and  (iii)  the  static  Coulomb  interaction  energy 

of  the  electrons  and  ions  in  their  oscillatory  equilibrium  positions. 

From  this  physical  model,  the  interaction  part  of  the  free  energy  is 

derived,  which  consists  of  a  quasi-lattice  energy  depending  on  the 

interaction  parameter  y,  and  the  free  energies  of  the  quanitized  electron 

and  ion  oscillations  (long  range  interactions).  Depending  on  the  degree 

of  ordering,  the  Madelung  "constant"  of  the  plasma  is  c t(a)  =  a  for  y>>l, 

Vo 

u(y)  =  ct  for  y  >  1,  and  a(y)  <ry  "  for  y  <<  1,  where  a  -  1  is  a  constant. 

The  free  energy  of  the  high-frequency  plasmonr.  (electron  oscillations)  is 
shown  to  be  very  small  for  y  >  1,  whereas  the  free  energy  of  the  low- 
frequency  plasmons  (ion  oscillations)  is  shown  to  be  significant  for 
y  >  1,  i.e.  for  proper  nonideal  conditions.  For  weakly  nonideal  plasmas, 
y  "  1,  both  the  electron  arid  ion  oscillations  contribute  to  the  free  energy. 
Unis,  novel  results  are  obtained  not  only  for  proper  nonideal  (y  >  1) 
but  also  for  weakly  nonideal  (y  <<  1)  plasmas.  From  the  general  formula 
lor  f  lie  free  interaction  energy  hr  of  the  plasma  for  0  <  r  ' ■  Y  ,  simple 
analytical  expressions  are  derived  for  AF  in  the  limiting  cases,  y  >>  1, 


1,  and  y  ■-<  1. 


INTRODUCTION 


In  the  classical  work  of  L)ebye  and  Htieckcl  on  electrolytes,  the  total 

Coulomb  interaction  energy  is  calculated  from  the  continuum  theoretical 

picture  of  every  ion  interacting  with  its  surrounding  space  charge  cloud. 

Using  more  sophisticated  methods,  similar  results  were  obtained  for 

1  2 

weakly  nonideal  plasmas  (y<<l)  by  Mayer  (cluster  expansion),  Ichikawa 

3  A 

(collective  variable  approach  ),  Vedenov  and  l.arkin  (graphical  density 
expansion),  and  Jackson^  (hydrodynamic  continuum  interaction  model). 

Based  on  different  methods  and  approximations,  investigations  of  moderately 

(y  >  1)  and  strongly  (y>>l)  nonideal  plasmas  were  given  by  Berlin 

6  7  8 

and  Montroll  ,  The  inter  and  Centry  ,  Kcker  and  Kroell  ,  Ebeling,  Hoffman 

9  in 

an  '  Kelbg  ,  and  Varobev,  Norman  and  Vilinov  ,  respectively. 

in  spite  of  differences  in  the  theoretical  approaches,  the.  leading 

terms  of  the  analytical  results  for  proper  nonideal  plasmas  (y>1) 

give  essentially  the  same1  formula  for  the  free  plasma  energy,  AF/NKT  =  -ay  + 

2  1/3 

blny  +  c,  due  to  Coulomb  interaction,  where  Y  =  Zc  n  /KT  is  the 

ratio  of  electron  -  ion  interaction  energy  and  thermal  energy,  and 

a,  b,  c  are  constants  depending  on  the  respective  approximations 

and  assumptions.  The  thermodynamic  functions  of  strongly  nonideal 

plasmas  (Y'-->l)were  also  determined  with  the  help  of  Monte  Carlo  and 

computer  methods  by  Brush,  Sahlin,  and  Teller^,  Hansen^,  I)eWitt^\ 

14 

and  Theimer  ,  respectively.  Although  computer  methods  provide  limited 
physical  insight,  they  are  useful  for  checking  the  quantitative  validity 
of  analvt ical  theories  . 


At  sufficiently  high  electron  densities,  for  which  Y  >  1,  classical 
statist  ical  theories  fail  due  to  thermodynamic  Instability^,  which  is 


inhibited  by  quantum  mechanics.  The  classical  plasma  pressure  would  collapse  for 
Y>  1  due  to  the  negative  electron-ion  interaction  energy,  whereas  in  reality 
the  pressure  remains  positive  in  a  plasma  due  to  the  Fermi  pressure  (exclusion 
principle)  of  the  electrons.  For  these  reasons,  we  present  herein  a  quantum- 
statistical  theory  for  nonideal  plasmas  based  on  concepts  similar  to  those 
used  by  Debye  for  solids^^.  The  application  of  this  model  to  proper 
nonideal  plasmas  (y  >  1)  is  justified  since  a  plasma  exhibits  a  quasi¬ 
crystalline  structure  for  y  >  0.1  before  it  undergoes  a  diffuse  transition 

into  a  solid,  metallic  state  at  a  critical  value  y  .  The  roll  of  the 

c 

longitudinal  phonos  of  the  Debye  theory  is  assumed  by  the  quanta  of  the 
plasma  oscillations  (plasmons)  in  the  case  of  the  quasi-crystalline  plasma. 

The  theory  is  also  applicable  to  weakly  nonideal  conditions,  since  the  quasi¬ 
lattice  energy  reduces  for  weak  ordering,  y  <■  1,  to  the  free  interaction 

energy  of  weaklv  nonideal  plasmas. 

The  theory  to  be  presented  takes  into  consideration  (i)  the  energy 
eigenvalues  of  the  random,  collective  electron  and  ion  oscillations  and 
(ii)  the  static  Coulomb  interaction  energy  (quasi-lattice  energy)  of  the 
electrons  and  ions  in  their  oscillatory  equilibrium  positions.  Thus, 
all  significant  long  and  short  range  Coulomb  interactions  are  considered. 

The  results  are  applicable  to  arbitrary  nonideal  plasmas,  0  <  y  <  Yc  ,  where  yc 
is  the  critical  ordering  parameter  at  which  a  phase  transition  into  a 
solid  metallic  state  occurs. 
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FtiYSlCAh  :o;:N|)ATIi>t;S 


Subject,  of  the  theoretical  considerations  are  quasi -homogeneous  high- 
pressure  plasmas  consisting,  of  electrons  of  charge  -e  and  density  n  =  N/V  and 

ions  of  charge  l/.e  and  density  u/Z  -  K/ ZV ,  with  typical  densities  in  the  range 

2()-3  IL (  _j  3  4  0 

IlT  cm  <  n  •_  lir  cm  and  temperatures  of  the  order-of-magnitude  T  10-  10  K. 

2 

For  these  conditions,  Lhe  Debye  radius  D  =  [4iinc  (1+Z)/KT]  11  is  D  =  6.901  x 

1  — B 

[ T/n ( 1+Z)  2  •;  10  cm,  i.e.  ,  is  smaller  than  the  atomic  dimension  and  the  number 
of  particles  in  the  Debye  sphere  would  be  N  4tniD  /3  <  1  for  Ik  10  cm  and 

>4  _3 

nv  10  cm  .  It  is  seen  that  the  concept  of  Debye  shielding  completely  breaks 

down,  and  statistical  theories  containing  the  Debye  length  as  a  characteristic 

parameter  would  be  physically  meaningless  for  high  density  plasmas. 

The  nonideal  behavior  of  plasmas  is  determined  by  the  interaction  parameter 

2  1/3 

'i  ,  which  is  the  ratio  of  the  Coulomb  interaction  energy  ~Ze  n  and  thermal 


energy  KT, 

y  =  Ze2n1/3/KT=  1.671  *  10“3Zn'/3/T.  <D 

It  follows  that  0.5Z  <  y  <  15Z  for  102^ctn  •_  n  <  102+cm  3  and  T  ~  10^  °K.  For 

Y  =  I,  the  nature  of  the  plasma  changes  from  a  "thermally  expanding"  (y^l)  to 
an  "electrostatically  contracting"  (y>l)  plasma.  For  y>l,  the  collapse  of  the 
plasma  due  to  Coulomb  attraction  between  electrons  and  ions  is  inhibited  by  Lhe 
Fermi  pressure  of  the  electrons,  i.e.  by  the  quantum  mechanical  exclusion  prin¬ 
ciple.  Thus,  in  the  region  0  <  y  <■  y  the  plasma  undergoes  a  diffuse  transi- 

e 

tion  from  a  nonideal  classical  plasma  (y  <  1)  to  a  quas i-crystall inc.  plasma 

(1  *■  y  ■  y  ),  with  an  incomplete  ordering  comparable  to  that  of  a  liquid, 
c 


An  understanding  of  strongly  nonideal  plasmas  has  been  attempted  via  the 

i  •  i  (6-14) 

model  ot  discrete  interacting  particles  in  a  dense  gas  -  .  lor  the  above 

reasons,  however,  it  appears  to  be  a  more  reasonable  procedure  to  calculate 


tlif  thermodynamic  functions  of  proper  non  idea  1  |>  J  a:  anas  i  i  om  Liu.-  picture  of 

col  Uu't  ive  1‘Ji‘i'lron  ;rwl  inn  oscillations.  In  this  app  roach ,  the  free  inter¬ 
action  energy  is  due  to  the  static  Coulomb  interaction  of  the  electrons  and 
ions  in  their  "equilibrium  positions"  (.‘ladeiung  energy)  and  their  oscillation 
energies  about  the  equilibrium  positions  (plasmon  energies). 

Since  the  plasma  volume  V  contains  ,\  electrons  and  N/k  ions,  there  exist 
SN  (high-frequency  branch)  and  iN/Z  (low  I  i  equeiicy  branch)  characteristic  fre¬ 
quent.  ies  'X .  ot  longitudinal  oscillations.  Kacli  plasma  oscillator  of  frequency 
U.  can  only  have  the  energy  (n.  +  '. )  h .  ,n  -0,1,2,...,  s>  that  the  energy  Fa  i  j 

l  tit 

of  a  plasma  state  with  n.  plasmons  ol  frequency  u .  is 

l  i 

K  {  i  }  “  )  n  .diui .  (2) 

I  i}  1  1 

where  {if  designates  the  entire  set  ot  given  eigeni nqaencies  m .  Ac¬ 
cordingly  ,  the  partition  function  Q  of  rlu-  longitudinal  plasma  oscillations  is 

t)  =  )’  e-K,i,/KT  =  II  l/d-e^i7^1').  (3) 

{ i  >  1U 


Front  Q,  the  thermodynamic  functions  such  as  the  pressure,  internal  energy, 
entropy,  etc.,  are  derived  in  Lite  usual  way,  e.g.,  the  free  energy  of  the 


plasmons  is 


F  -  -KT  In  t)  =  K.T  ^  ln(l-e  , 


In  the  limit  V*-"  ,  the  discrete  eigenf  reipienc  ies  o> .  are  t'ep  Laced  by  con¬ 
tinuous  ones,  ti>=id(k)  ,  in  accordance  with  the  uispersion  law  for  space  charge 
waves  of  wave  length  =  2 it / k ,  ()  y  k  <  k. 

I.  Fleet ron  Oscl  llations.  The  high-frequency  branch  of  the  space  charge 
waves  i  s  cltio  to  longitudinal  electron  oscillations.  Their  i  requency  m  is  for 
classical  (n  '•  n)  and  completely  degenerate  (n  •  n)  electrons  given  by-*-^ 


ur  [1  f(r  /4tt)ZY  '  (kr  )  ) 
]>  e  v 


1 


n  •  *  n 


(6) 


<7 


or  -  m 


11  +  zsr  (-^)1/|<t>7lzi'‘lik':u)  i. 


where 


.  y, 

n  =  2(2)imKT/lr  ) 

1/ 

,,i  =  (Anne-  /m)  ‘ 

P 

r-  =n-^ 

t? 


(7) 

(H) 

(9) 


are  the  critical  electron  density,  the  plasma  frequency,  and 

the  mean  electron  distance  (k  =  c  /c  of  the  electrons,  and  m  is 

e  p  v 

their  mass).  Since  k  ~2ir/r  (oscillations  with  \  <  r  are 

max  e  e 

physically  inconceivable),  the  electron  oscillations  propagate, 
w  =  w(k)  >  w  ,  in  non  ideal  plasmas. 

2.  Ion  Oscillations.  The  low-f requency  branch  of  the  space 
charge  waves  is  essentially  due  to  ion  sound  waves.  Since  the  ions 

are  presumed  to  be  nondegenerate,  the  frequency  of  the  ion  oscillations 

,  17) 
is  given  by-- 


1/ 

u  -  A  (k) (k  .KT/K)  ‘  k 


(10) 


where 


'  (k) 


1  + 


z(kq/<  ; ) 


VA 


1  +  (k:  /Air)ZY  '  (kr  )' 


(11) 


<5  ( k )  =  1  ,  n  •  n  ,  (1.2) 

is  a  correction  factor  of  magni tue-of-order  1,  which  shows  the 

influence  of  the  electrons  on  the  ion  oscillations  (M  =  mass,  r.  = 

c  /c  of  the  ions). 

P  v 

In  weakly  nonideal  plasmas,  y  •  i,  the  electron  sound  waves  are 


strongly  damped  for  wave  lengths  X  •  D,  due  to  trapping  of  the  resonance 


electrons  with  thermal  speeds  comparable  to  the  wave  speed.  For  prope 
nonideal  plasmas,  y  %  1,  the  number  ol  particles  in  the  Dobvc  splure 

,  -H 

4nDJ/3  is  no  longer  large  compared  with  one  and  D  ■■  10  cm  is  smaller 
than  atomic  size,  so  that  thermal  Landau  damping  is  no  longer  feasable 
For  this  reason,  electron  oscillations  should  exist  for  wave  Length 


VI 


ST  AT  I  ST  I  f’.AI .  TIIIT.MUDYN AM  I  (  S 

In  the  plasma  iimli'r  cons  Ldernt  ion,  I  lie  e  i  oil  run:;  aiui  ions  interact 

through  their  longitudinal  Coulomb  fields  (transverse  electromagnetic 

2 

interactions  are  negligible  for  KT  me  ) .  1 1  u  •  electrons  (s  -  e)  and 

ions  (s  =  i)  have  tliermal  velocities  c  and  random  collective  mean 

mass  velocities  v  due  to  their  oscillatory  wave  motions  about  the 

s 

equilibrium  positions,  so  that  /;/  i  f  d "V  =  f),  and  •  v  -■  =  0  where 

I  is  the  equilibrium  velocity  distribution  of  the-  species  s  and 
s 

designates  a  spatial  or  temporal  average.  The  resulting,  Hamilton  function 
with  Coulomb  interaction  leads  to  a  tree  energy  of  the  plasma  of  the 
form : 

1’  )'  h'J0)  +  iiM  +  l !s  .  (13) 

s=e , i  s-o ,  i 


1^°^  is  the  ideal  tree  energy  of  tin  non  im  eraet.  i  ng  plasma  components 
s.  K  is  the  Coulomb  interaction  energy  ol  the  electrons  and  ion';  in  their 
equi librium  positions.  F^  .  is  the  free  energy  of  the  electron  and  ion 
oscillations,  i.e.  of  the  high  and  low  frequency  plasmons,  llq.  (A). 

It.  should  he  noted  that  bq.  (’>)  takes  into,  eons  iderat  ion  all 
significant  short-range  and  long-range  Coulomb  interactions  by  means  ol 
the  Made  lung  energy  F.,,  and  the  plasmon  energies  F  .  As  is  evident  from 


17) 


the  derivation  of  bqs.  (5)-(6)  and  flO),  in  wiiieh  terms  of  order  m/M  are 

neglected  compared  to  1,  Fq .  ((1)  coin  ains  the  e-e,  e-i,  and  i.-i  Coulomb 

-1/1 

interactions  at  distances  A  >  n 

1.  Free  Energy  Fs°^ .  In  high  pressure  plasmas,  the  electrons  are 

153/2  -3 

part  iallv  degenerate  for  dens i tie:  a  where  n  =  A. 828  x  10  T  | cm  ], 


whereas  the  ions  behave  in  general  classically. 


for  the  free  energy  of  the  ideal  electron  gas 


8) 


Fermi  statistics  gives 


4<l 


i'^'0  =  -  NKT  U.1/.,(p/KT)/U1/  (n/kT 


(14) 


whei\ 


U  (i./KT) 


J _  C  _ x  P  dx _ 

(i’  +  l)  J  ^x-m/KT  f  ! 


I/:’.  •*/:? 


(13) 


and 


n  =  2  (2irmK T/h“)  32“  U^.,  (p/KT) 

19) 


(16) 


defines  the  Sommer  f  eld  integrals  ,  .and  let  <-rin  i  lies  tie*  eliemicnl 

potential  p=p(n,  T)  of  Lite  eleetrons,  respectively.  The  free 

energy  of  the  translational  degrees  of  freedom  of  I  lie  classical, 
18) 


ideal  ion  gas  is 

F(o)  =  -  (N/Z)  KT  ln[ (2nMKT/h2) 3/2 
\ 


(17) 


2.  Quasi-Lattice  Rnergy  E^.  The  equilibrium  positions  of  the 

electrons  and  ions,  about  which  the  electrostatic  oscillations 

occur,  form  an  electron  "lattice"  and  an  ion  "lattice",  with 

18) 

an  incomplete  ordering.  By  means  of  the  Ewald  method  ,  one 
calculates  the  Coulomb  ineract ion  energy  of  the  quns i-cubic al ly 
centerid  elect  run- ion  lat  tices  as 

Km  =  -  re,  NKT  ,  <»  =  ct  =  1.451  for  Y  -•  l.  (18) 

As  the  ordering  of  the  plasma  increases  with  y,  a  ( y )  is  a  weak 

function  of  y  such  that  asymptotically  rx  =  1.541  for  y  1.  lr.q.(18) 

2  _ 

indicates  that  -E  /N  Ze  / r  Is  of  the  order  of  the  average 
M 

e- i  interaction  energy.  For  weak  ordering,  y  ■  1,  it  will  he  sho.Ti 

V, 

Y 


that  <t 
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3.  High-Frequency  Contribution  F  .  Since  the  number  of  Longitudinal 

modes  with  wave  numbers  between  k  and  k  +  dk  in  volume  V  is  V4tt  k2  dk/ (2tt  )  2  , 

Eq.  (4)  gives  for  the  free  energy  F  of  the  high-frequency  electron 

oscillations  of  energy  Tim (k) 

ka 

Fo/KT(V/2tt2)  =  l!n{  .l-cxpl-ti'to(k)/KT)  }  k:'dk  (19) 

0 


where 


,  ,  v, 

m(k)  =  u  (1+a2  k2 )  '<  , 


cz/ur  =  (k  /47T)(Z/>)rz  ,  n  <<  n  , 
m  p  e  e 

(3/5)  vJf/u2  =  (-~)  ^(-v-)  (Z/y)V  ,  n  >>  n  , 

1  p  2Uti  on  e 


(20) 

(21) 

(22) 


by  Eqs.  (5) -(6).  The  speed  of  sound  and  the  Fermi  speed  Vp 


of  the  electrons  are 

c  =  (k  KT/m) 1/2  , 
m  e 


v,,  =  'ff(3Tr9n)  ^/m. 

r 


The  number  of  modes  in  (0,  kfi)  and  V  equals  the  number  3N  of 
degrees  of  freedom  ot  the  electron  gas,  i.e.. 


(2r)  ’V 


K  1/ 
rc  4irk2dk  =  3N,  k  =  (18irn)  h 

J  e 

0 


(23) 


(24) 


Integration  of  Eq.  (19)  by  parts  yields,  under  consideration  of 
k!  KTV/hn2  =  3NKT,  for  the  free  energy  of  the  high-frequency  plasmons: 


Tim 


3NKT  bill  -  exp [  - 


JL 


KT 


(  1+a'  k^  )  ^  |  f  -  F  ( 


Tim 


KT 


■  •lk0) 


(25) 


who  re 


lito  hi.',  ak  i(/,  ,  ;■ ,  -  ’/'  , 

ak  >  ■  — ' c<ak  J  "  - 

KT  KT  (tfco  /KT)  (i+x2)  1 

0  ° 


(261 


and 


li'.i  /KT  =  ( 4  ;t )  ^  (A  /r  )  ( y / i . )  ^ 

p  e  e 


>  =  b/(mKT),/?  , 

o 


(27) 
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ak  =  k  ^  (9it  ^?/4)  ^(Z/y)  ^  ,  n  <•-  n  , 
0  0 


ak£  =  2  ^&7i  ( 3 / 2 )  (n/n)  ^(Z/y)  ^  ,  n  n. 

By  means  of  the  successive  substitutions,  (f)  x-sinln,,  dx  =  cosh  c,  df 
and  (ii)  c  =  ("Rio^/KT)  coshf,,  dc  »  (Roj^/K'1)  slnh  fd(.,the  integral  (26) 
is  transformed  to 


F(ep’ake)  =  (ake':p)-3  \  (f2_t?p)  (c'  -l)  ‘  <l’ 


Vo  /  _ 


where 

e  =  diw  /KT,  e  =  e  cosh(arcsinh  ak  ).  (3! 

p  p  ’  e  p  e 

Since  the  leading  expression  in  Eq.  (25)  is  the  logarithmic  term,  it 

is  sufficient  to  give  for  F(e^,  a^e)  tl.e  series  approximation  (Appendix), 


,3.  -%  , 


F(c  , ak  )/2  '2  (ak  )  3e 
P  e  e  p 


l  e  mEp  l  („)(2e  )  n  m  +  n)y(^  +  n,  (c  -e  )m),  e  <  3e  ,  (32) 

,=  i  n=n  n  P  2  e  p  e  p 


where 


5  e  -e  3 

,  „  -  +  n  eP  2+n_ 

Y (-  +  n,  (e  -e  )ra)  =  m  /  u  e  mudu 

2  e  p  J 


is  the  incomplete  gamma  function —  .  Since  in  general  y/Z  ~  1  for 


Cp  •-  e^  <  ^Ep,  the  expansion  (32)  is  useful  where  simple  approximate 
relations  do  not  exist. 


43 


4.  Low  Frequency  Contribution  F  With  the  number  of  modes  in  the 
interval  dk  at  k  and  volume  V  given  by  V4rk2dk/ (7  it ) 3 ,  Eq .  (4)  yields  for  the  free 
energy  of  the  low-frequency  ion  oscillations  of  energy  Tfm(k) 


F./KT(V/2tt2) 


In  { 1-exp  [  --Rw  (k )  /KT  ]  >  k2dk 


where 

u>(k)  =  6(k)cMk, 

1/ 

cM  =  (r.KT/M)  /?  . 

by  Fqs.  (10)  and  (12).  The  number  of  inodes  in  (0,  k^)  and  V  equals 

the  number  3N/Z  of  degrees  of  freedom  of  the  ion  gas,  i.e., 
k 

(2n)~3V  (  4Tik?dk  =  3N/Z,  kf  =  (18u?n/Z)  ^ 


Partial  integration  of  Eq.  (34)  gives,  under  consideration  of 
k?  KTV/6ir2  =  3(N/Z)KT,  for  the  free  energy  of  the  low-frequency  plasmons: 


F.  =  3 (N/Z)KT  I  fn {1-exp f- 


6 (k . ) k . ] )  -  G(k. 

i  t  j 


where 


(l(k  )  =  k  -3  f  1  U(k)4k6‘  (k)lk3dk 

1  KT  1  j  e(fiCjij(KT)6(k)k 


Since  the  dispersion  factor  <5(k)  is  a  bounded  function  varying  very 

1 L 

little  with  k  such  that  1<  S(k)  ^  (1+Z)  '  for  k  r (0,  k^),  6(k)  can 

be  approximated  by  an  average  value  6  , 

6  (k)  =  6  1,  n  ^  n  . 


Since  in  addition  the  logarithmic  expression  is  the  dominant  term  in  Eq.  (38), 
the  integral  (39)  is  approximated  by 


G(e  )  =  e  f1  t'(ef  -1)~ 1  dr 
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where 


e  =  "hc^Sk/KT,  e  L  * -fic^k  /KT. 


70) 


G(e^)  has  the  serai-convergent  series  expansions, 

G(Ci)  =  3U— +  20zi  +-  *  -1  ’  ei  <<  lf 
G(6l)  =  -jy  +  0[e  c i] ,  ei  >>  1. 

This  completes  the  formal  mathematical  aspects  of  the  theory, 
the  physical  implications  of  which  require  further  elaboration. 


(42) 

(43) 

(44) 
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APPLICATIONS 

For  applications  of  the  theory  to  strongly,  Intermediate,  and  weakly 

nonideal  plasmas,  it  should  be  noted  that  the  dimensionless  parameters 

y/ Z , 'flu)  /KT,  ak  ,  and  n/n  occuring  in  Eq.  (25)  for  the  free  energy  F 
p  e  e 

of  the  high-frequency  plasmons  can  not  be  varied  independently.  Since 
both  y/Z  and  A__/r^  increase  with  increasing  n  and  decreasing  T, 

Tho  /KT  ~  (Ag/r g)  (y/7.)  '2  varies  over  a  large  n-T  region  similar  to  (y/Z)1, 

Eq.  (7).  Numerically, 

y/7  =  1.670x  10“3n1/3/T,  -fito  /KT  =  4 . 310x1 0_7n ^  /T ,  n/n  =  2 . 071xl0”16nT"^ 
ak>  =  2 . 219  (y  /Z)  12  ,  n  <<  n 

akg  =  1. 910(n/n)  73  (y /Z)  %  n>>ii.  (45) 

E.  g.,  for  T=104°K,  y/Z£  1  if  n  <’  1021cm“3  and  !•:,»/ KT  >  1  if 
n  >  5xl02°cm~3.  For  T=103°K,  y/Z  <  1  if  n  i  10KVm"3,  etc.  Thus, 
for  typical  conditions  of  nonideal  plasmas  y/Z  and  Urn  /KT  are  of  the 
same  order  of  magnitude.  It  is  also  recognized  that  in  general 
n/n  •■>  1  if  y/Z  >>  1,  and  n/n  <<  1  if  y/Z<<  1. 


Tn  Eq.  (38)  for  the  free  energy  F^  of  the  low  frequency  plasmons, 
only  one  characteristic  parameter  e^  occurs  since  6(k)~  4  ~  1.  By  Eq. 
(42) ,  this  parameter  is 


(18ir2)1/3Ki1/2  5  - —  =  2.158x10  5Z  ^  (•£}•)  ^  ^  ^  K<  1  (46> 

r . 

i 


where 


= -ft/ (MKT) ^  ,  vi  =  (n/Z)^3  .  (47) 

Accordingly,  for  typical  nonideal  plasma  conditions,  it  is  <<  1 
since  >  /r^*-'  1  (classical  ions)  although  in  general  A^/r^  > 
electrons)  for  y/7,  >  1  or  1i«)  /KT  >  1. 


1  (degenerate 


46 


The  deviation  AF  of  the  free  energy  of  a  nonideal  plasma  from 
ideality  is  by  Eq.  (13)  due  to  the  quasi-lattice  energy  F.^  and  the  plasmon 
energies  . , 

AF  =  Ew  +  Y  F  .  (48) 

s=e ,  l 

Since  the  theory  of  electron  oscillations^^ has  not  yet  been 
developed  for  arbitrary  degrees  of  degeneracy  (n  >  n),  the  contributions 
of  the  electron  iscillations  to  AF  in  the  cases  n  n  and  n  ?  n  have  to 
be  estimated  from  the  dispersion  equations  for  n  <<  n  [F,q.  (5)]  and 
n  >>  n  [F.q.  (6)],  respectively.  Fortunately,  it  turns  out  that 
]F  | < < | AF |  for  y/z>l>  so  that  quantitatively  rolyable  approximations 
for  AF  can  be  derived. 


1.  Strongly  Nonideal  Plasmas.  By  Eq.  (6)  the  spectrum  w(k)  of 

electron  oscillations  extends  over  a  band  Am  -  w  above  the  plasma 

-  <  -  P  %  -3 

frequency  for  y/Z  »  1  since  krg-  k^r  ~  1  and  (n/n)  ~  1.  Application 

of  the  mean  value  theorem  for  integrals  to  Eq .  (25)  shows  that  the  free 


energy  F^  of  the  high-frequency  plasmons  vanishes  exnonentially  for 
£  i .  e .  y  /  7.  -*  oj  : 


?  /3NKT  =  (ln{ i  -  exp [ -c  (l+a2k  2 )  ^  ] 
e  y  p  e 


e  (ak  ) 
p  e 

exp[Cp(l+x‘ 


— 77 -  r  x^l+x2)  ^dx^  v  0,  c 

),/2]-l  J  /  p 


0  -  x  -  ak 

e 


Accordingly,  |  |  / 3NKT  < <  1  for  £^>>1,  i.e.,  y/7,  •>!.  On  the  other 

liand,  the  free  energy  of  the  low  frequency  plasmons  is  by  Eq.  (38)  for 


nondegcuerate  ions 
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F.  =  3(N/Z)KT[£he,  -(1/3)]  = 

^  ^  V 

3(N/Z)KTU'hy  +  £n[  (18tt2/Z4)1/3  ?  ]  -  (1/3)},  e.  <■<•'  1.  (50) 

e?/fi  1 

It  is  noted  that  y/Z  >>  1  is  compatible  with  e.  =  1Tc^k^6/KT<<  1  as 
explained  above. 

Equations  (49)  and  (50)  demonstrate  that  the  contribution  of  the 
electron  oscillations  to  the  free  energy  is  negligible  in  strongly 
nonideal  plasmas,  y/Z  >>  1.  In  this  limit,  the  nonideal  part  of  the 
free  energy  is  due  to  the  quasi-lattice  energy  E ^  and  the  ion 
oscillations , 

AF/NKT  =  -  ay  +  (3/Z)6ly  +  (3/Z)£rt  (f^A^)  -  0/7.),  y/Z  »  1,  (51) 

where 

v B  =  e2/h,  6  =  (1  8tt2 Z~^)  .  (52) 

Note  that  tny  depends  on  both  n  and  T  whereas  fiiPc^/Vg  depends  only 

Q  1  I 

on  T,  where  the  Bohr  speed  is  vfi  =  2.118  x  10  cm/sec  >>  cM  =  (i^KT/M)  '2 . 
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It  is  remarkable  that  the  electron  oscillations  contribute  little 

to  the  free  energy  compared  to  the  ion  oscillations  for  y/Z  •>  1.  This 

result  holds  even  for  moderately  nonideal  conditions,  y/Z  >  1.  Thus, 

we  disagree  with  the  formula  "F  =  ne  +  3NKTf’lt  (1'uj  /KT>"  stated  without 

o  o 

21 

derivation  for  nonideal  plasmas  by  Norman  and  Starostin-  ,  according 
to  whom  "all  the  vibrations  have  exactly  the  same  frequency  u 
near  the  plasma  frequency  <0^".  The  derivation  of  this  formula  requires 
fiw(k)/KT  <<  1  for  the  electron  oscillations,  which  implies  y/Z  <<  1, 
but  the  latter  inequality  contradicts  their  assumption  w(k)  =  wo  =  w^, 

since  the  frequency  spectrum  extends  over  a  large  band  Aw  >  w  above 
0)^  for  y/Z  <<  1.  For  these  reasons,  the  free  energy  proposed  by  them 
is  not  applicable  to  proper  nonideal  plasmas,  y/Z  '  1,  nor  is  it  correct 
for  weakly  nonideal  conditions,  y/Z  <<  1. 

2.  Intermediate  Nonideal  Plasmas.  For  intermediate  nonideal 

conditions,  1  .  y/Z  <  10,  the  spectrum  m(k)  of  electron  oscillations 

extends  over  a  region  Aw  0[w  ]  above  w 

P  P 

k r  1  k  r  ~  1.  Also  in  this  case,  a  relatively  simple  formula  can  be 

e  e  t 

devised  for  the  free  energy.  The  logarithmic  term  in  F  ,  Eq .  (25) 

is  negligible  compared  to  that  in  F  ,  Eq.  (38).  for  y/z '■  1  since 

t  .-.Tic  iSk./KT  for  y/Z  >  1  bv  Eos.  (45)  and  (46),  respectively, 
p  Mi 

Accordingly,  the  nonideal  part  (48)  of  the  free  energy  is  for  inter¬ 
mediate  non  ideal  plasmas: 

E/NKT  =  -7ty  +  (3/Z)fuy  +  ( 3/Z)  f«  (8.^/v^  -  (3/Z)0(r 
-  3f:  (t  ,ak(j )  ,  y/z  T  1. 


by  Eq.  (6) 


(n/n) 


(53) 
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For  )/'/.  -  1,  the  ions  can  bo  assumed  to  bo  non-do^cno  rato  ,  t  ^  -  fic^ok  /KT  ■  1 


by  Eq.  (4b),  so  that  the  ion  integral  (41)  rodnoos  to 

«(>'.)  =  i/i.  r .  1  . 

Since  t  ^  r  1  and  ak^t  r.  1  [  Kq .  (Ar,)]  for  y/Z  '  1  >  ti,t.  electron 
integral  (U))  is  significantly  smaller  than  l;(r  .)  =  1/3, 

0  •  I  (,  ,ak  )•  (z2~z?)^2  (ak  e  )  ~  J  f  f  f  ( I-c*"  e°/ 1-e P  °)  --  1,  y/Y.i  1. 
peep  e  p  1 

The  lower  and  upper  bounds  of  F(i_^,ak  )  have  boon  obtained  by  means  of 

the  mean  value  theorem  for  the  integral  (30), 


( (c  ,ak  )  =  (ak  e  )' 
p  e  e  p 


V  ^ 


-l)d> 


(14) 


(05) 


(56) 


While  for  strongly  nonideal  conditions,  the  contribution  of  the 

electron  oscillations  to  the  free  energy  is  completely  negligible, 

this  contribution  is  still  insignificant  for  intermediate  nonideal 

conditions,  y/Z  >  1,  by  Eq.  (55).  For  more  exact  evaluations,  the 

small  term  F(c  ,ak  )  in  Eq.  (53)  can  be  computed  from  Eq .  (30)  or  (32). 
p  e 


3.  Weakley  Non  ideal  Plasmas.  Although  the  theory  of  weakly 

1-5) 

non  ideal  systems  is  well  understood,  it  is  interesting  to  investigate 
whe i her  the  present  model  for  proper  nonideal  plasmas  gives  reasonable 
results  in  the  limit  y/Z  <<  1.  For  y/Z  <<  1  it  is  ak^  >>  1  by  Eq .  (45), 
and  the  spectrum  w(k)  of  electron  oscillations  extends  over  a  large 
region  Am  '•>  u>  above  w  by  Eq.  (5).  The  electron  integral  becomes  for 


ak  >  >  1 , 


F(e  ,ak  )  =  c 

P  e  P 


a  k 

(ake)"3J  (eCPX-l)_1x3dx,  y/Z  « 


1, 


(57) 


r(,  ako)='/i|,-l(1  ak^)  F  i(r  akp)2-  ...],  r  .*  «  1.  (58) 
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.Ait itougli  t  ak^  is  independent  of  y/Z  by  Eqs.  (27)  and  (23),  the  expansion 
(58)  is  valid  since  the  electrons  are  certainly  nondegenerate,  }■  /  r  <"•'  ] 


for  y  /  7.  <<  1 ,  and 

1/  1/  1/ 

e  ak  =  (  4  )  •  (9n  2/4)  3  X  /r  1,  >  /  r  <•  1. 

p  e  e  e  e  e  e 


(59) 


For  nondegenerate  ions,  the  integral  (41)  is  C(t^)  =  1/3  by  Eq.  (43) 
since  e.  --v  1.  Thus,  one  obtains  from  Eqs.  (18),  (25)  and  (38)  for 
the  interaction  part  of  the  free  energy  ol  weak  1 v  nonideal  plasmas: 


AF/NKT  =  -u(y)y  +  (3/Z)i';ty  +  (3//.)  tn  (dc  /v  ) 

+  36l(epake)  -  (1+Z_1),  y/7."  1, 

where  t. bo  logarithmic  term  in  Eq.  (25)  has  been  o-  ponded  for 


(60) 


c  ak  ■ *  1 . 

P  e 

In  Eq.  (60),  u(y)  is  the  Madelnng  constant  of  the  weakly  nonideal 
plasma  with  weak  electron  and  ion  ordering,  n(y)  *0  for  y  -»  0.  Comparison 
of  the  term  -a(y)y(NKT)  in  Eq.  (60)  with  AF  = 

-  (NET)  (2/3)  it  ^  (1+Z)  ^e  3n  '  ‘  (KT)  ^  of  the  Debv e-Hueckel  theory—^ 

(weakly  nonideal  plasmas)  yields  the  result 

<i(y)  =  (2/3)ir!/>  (l+Z~V/?  yl/2  ,  y/7.  <  1  .  (61) 


The  previous  theories  of  weakly  nonidoal  plasmas  do  not  lead  to  the 
logarithmic  terms  in  Kq .  (60)  since  they  ’o  not  take  into  account  the 

effects  of  electron  and  ion  oscillations. 

The  presented  theory  is  applicable  to  nonideal  plasmas  in  the  gaseous 
phase  (0-y-l)  and  the  quasi-liquid  (l-v-y,)  phase.  Whether  the  Eqs.  (25) 
and  (38)  are  applicable  to  hot  (T>10’CK)  plasmas  in  the  solid  phase 
( y  y  )  cannot  be  judged  at  this  time,  since  not  enough  is  experimentally 
1  own  about  the  latter,  extreme  state  ol  matter. 


APPENDIX:  Expansion  of  I  (t.  ,ak  ) 

1>  <• 

The  integral  (23)  is  conveniently  rewritten  in  the  fori: 

f(e  ,ak  )  =  (ak  e  )  I  (e  ,c) 
p  e  e  p  p 


where 


lV‘>  j  <  •>-e)%(ee1)-;J<  ,  o-  ,.p 


Since  r  >  0,  i.e.  e  <  1,  there  exists  the  series  expansion, 

(e* -1 )  1  =  £  e  mC  ,  e  >  0  . 

m=l 


The  snhst  itntion ,  u  = c  -  e  .  du  =  dt  ,  and  Kq .  (A3)  tiansiorm  Kq .  (A2)  to 


Kr, ,r)  =  l  e  m  P 


me  r  V  h  -mu 


r  •/ 

I  u 


(u+2ep)  (•  du. 


For  u <  2e  ,  i.e.,  e  <  3e  ,  the  binomial  expansion, 

P  P 

(u+2ep)^  -  (2cp)^  l  u/2cp<l,  (/ 

n=0 

is  used,  which  reduces  Eq.  (A4)  to  the  double  series: 

I  (t-  ,e)  =  (2l  )  l  e“m'p  l  (  n  )  (2r  )~n  m  _(7 +n)Y(4- +  n,  (c-e  )m)  , 


m=l  n~0 


c  <  3r  , 

P 


where 


5  .  e-e  3  . 

-z-  +  n  .  p  -rr+n 

5  ,  ,  ,  2  f  2  -|nu  , 

V  (— +n  ,  (e  -c  )m)  =  m  i  u  e  du , 


is  the  incomplete  gamma  function,  which  is  tabulated.  In  an 


amologous  way,  the  integral  (A2)  can  be  solved  for  u  >  2ep,  i.e., 
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V.  Variational  principles  and  canonical  theory 

FOR  MANY  -  COMPONENT  PLASMAS  WITH  SELF-CONSISTENT  FIELD 


ABSTRACT 

The  Lagrangian  L  «*  ///  L( q  ,  Dq  ,/D t,  Y*q  )d3r  and  Hamiltonian 

8  S  8 

H  ///  N(p  ,  q  ,  V-q  )d3r  for  many-component  plasmas  with  self-consistent  field 
s  s  s 

interactions  are  derived.  The  canonical  fields  p  (r,t)  and  q  (r,t)  are  defined 

S  3 

by  p  =  m  v  and  Sq  /St  =  n  v  ,  respectively,  where  n  (r,t)  is  the  density  and 

S  8  S  3  S  8  S 

■  ►  — ► 

v  (r,t)  is  the  velocity  of  the  component  s,  s  *  1,2,... N.  Based  on  the  action 
principle,  the  Lagrange  and  Hamilton  equations  of  motion  for  the  components  (s) 
are  presented  as  functional  derivatives  of  L  and  H  with  respect  to  the  canonical 
momenta  p  and  q^.  it  is  shown  that  the  new  formulations  oi  many- component 
plasma  dynamics  have  mathematical  advantages  compared  with  the  conventional 
variational  principles  and  hydrodynamic  equations. 


54 


INTRODUCTION 

Since  the  work  of  Clebsch  —  it  Is  known  that  the  space  and  time  dependent 
dynamics  of  perfect  fluids  can  be  derived  from  variational  principles.  - 
To-date,  iiowever ,  no  variational  principle  has  been  found  the  extremals  of 
which  give  the  Navier-Stokes  equations  with  viscous  dissipation.  A  comprehensive 
review  of  the  most  significant  contributions  to  variational  principles  in  fluid 
mechanics  has  been  presented  by  Serrin.  —  Similar  variational  principi  s  for 

perfect  one-fluid  plasmas  and  magnefohydrodynamic  fluids  have  been  given  by 

7  8 

Newcomb  —  and  Zakharov.  — 

Herein,  a  new  variational  principle  and  canonical  theory  for  perfect,  non- 

relativistic  many-component  plasmas  is  developed,  in  which  the  charged  particles 

interact  through  the  self-consistent  electric  field.  Interactions  through 

.••elf-consistent  magnetic  fields  are  neglected  for  non-relativistic  velocities 

v  i.  and  temperatures  kT  <<  m  c2  of  the  electrons.  We  show  that  the  electro- 
e  e  e 

hydrodynamic  equations  of  the  components  s  =  L,  2,  ...  N  of  the  plasma  are 

9 

derivable  as  functional  derivatives  —  of  the  Lagrange  function 

q  -V  ,  --V  41 

L  =  ///  L(r,t)dJr  and  Hamilton  function  H  =  Iff  n(r,t)d  r  of  the  plasma  of 

it 

volume  il.  The  Lagrange  L  and  Hamilton  tt  densities  are  functionals  of  the 

"►  ">■  . 

canonical  congugate  variables  and  defined  by  =  m^v^(r,t)  and 

8q  /Dt  =  n  (r,t)v  (r,t),  where  m  is  the  mass,  n  is  the  density,  and  v  is  the 

S  S  S  S  S  s 

(local)  average  velocity  of  the  particles  of  the  s-th  component. 

The  variational  principle  and  canonical  theory  for  many-component  plasmas 
provides  a  theoretical  foundat Ion  for  the  analysis  of  classical  phenomena  such 
as  nonlinear  electron  and  ion  acoustic  waves,  space  charge  fluctuations, 
electrostatic  turbulence,  and  e lect rohydrodyn.imi c  init ia l-boundary-vn 1 ue 


problems.  The  Hamilton  function  and  canonical  formalism  were  derived  as  a 
theoretical  basis  for  a  statistical  mechanics  of  hydrodynamic  plasma  fields 
concerned  with  the  evaluation  of  the  free  energy  of  electron-ion  plasmas,  which 
are  nonideal  due  to  the  sel  f-consist  ont  field  interactions.  The  statistics  of 
fields  and  wave  modes  represents  an  alternative  to  the  treatment  of  many-body 
systems  within  the  framework  of  the  discrete  particle  picture. 
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VLASOV  PLASMA 


The  subject  of  consideration  is  a  plasma  consisting  of  an  electron  component 

and  N-l  ion  components,  which  interact  through  the  self-consistent  electric  field 

E(r,t)  of  all  charged  particles.  The  hydrodynamic  equations  for  each  component 

s  of  the  plasma  are  obtained  as  moments  of  the  Vlasov  equation  for  its  velocity 

distribution  function  f  (v,r,t),  and  are  in  absence  of  collisions  (s=l  ,2  , .  . .  N) 

s 

7—  (n  m  v  )  +  V ■  (n  m  v  v  )  =  -  VP  +  n  e  E,  (1) 

3t  s  s  s  ssss  s  s  s  ’ 

|—  n  +  V- (n  v  )  =  0  ,  (2) 

3 1  s  s  s 

|-  P  +  V •  (P  v  )  =  -(->  -1)P  V-v  ,  (3) 

Sts  ss  s  ss 


where 


V-E  =  e"1  V  n  e  (4) 

o  u ,  s  s 
s=l 

->  ■>  ->■ 

relates  the  self-consistent  field  E(r,t)  to  its  space  charge  sources  n  (r,t)e  . 

s  s 

(c  is  the  dielectric  permittivity).  The  charge  of  the  particles  is  e  =  -e  for  the 
o  s 

electrons  and  =  Z^e,  |z  |  =  1,2,3,...,  for  the  ions  (e  is  the  elementary  charge), 
The  moments  are  the  density,  velocity,  and  pressure  fields  cf  the  s-th  component. 


ng(r,t)  =  ///  fg(v,r,t)d3v  , 

— <X) 

-foo 

v  (r,t)  =  n  1  ///  vf  (v,r,t)d3v  , 

s  s  s 

—00 

+°° 

P  (r,t)  =  (2/3)  III  %m  (v  -v  )2f  (v,r,t)  d3v. 


P,  v  ... 


Furthermore,  m  is  the  mass  of  the  s  particles  and  =  c  /c  is  the 

S  &  So 

adiabatic  coefficient  of  the  s  component  (v  is  the  particle  velocity). 

Equations  (1)  -  (4)  describe  the  macroscopic  dynamics  of  an  N-component 
Vlasov  plasma,  —  in  which  external  magnetic  fields  are  assumed  to  be  absent, 
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and  internal  magnetic  fields  and  binary  correlations  are  negligible. 
Maxwell's  equations  indicate  that  in  the  absence  of  magnetic  fields,  H  =  (5: 


V  x  E  =  0 


-r 

;»E  v  -*  A 

—  +  )  n  e  V  =0 

It  .  S  S  S 


Thus,  the  self-consistent  field  K(r,t)  is  longitudinal,  and  the  displacement 

(e  3E/3t)  and  convection  (n  e  v  )  currents  annihilate  themselves  at  every  point 
o  s  s  s 

(r,t).  This  effect  was  first  noted  by  Ehrenfest  in  his  discussion  of 
"accelerated  charge  motions  without  electromagnetic  radiation".  —  Typical 
examples  for  accelerated  charge  motions  without  transverse  radiation  field 
are  electron  and  ion  acoustic  waves,  and  space  charge  turbulence. 

Substitution  of  Eq .  (b)  into  Eq.  (4)  yields  the 
conservation  equation,  0  £  ^e^/lt  -f  V-^ne^v  =  0,  fcr  the  space  charge.  The 
latter  equation  is,  therefore,  not  independent  [in  addition,  it  is  obtainable 
by  summing  Eq.  (2)  multiplied  by  with  respect  to  s) . 

Due  to  the  absence  of  i)  external  magnetic  fields  and  ii)binnry  interactions 
(vanishing  viscous  and  thermal  dissipation) ,  the  pressure  in  Fqs.  (1)  and  (3)  is 

isotropic.  The  closure  of  the  Eqs.  (l)-(4)  is  accomplished  by  setting  all 

higher  order  moments  of  the  velocity  distribution  zero,  such  as  the  off- 

diagonal  stress  components,  n„  ~  0,  the  third  order  heal  flux  tensor, 

q.  -  0,  etc.  The  resulting  electrohydrodynamic  or  hydrodynamic  Vlasov 
t  J  k 

equations  are  generally  accepted  in  the  literature. 
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EXTREMUM  PRINCIPLE 


The  dynamics  of  the  many-component  plasma  is  now  condensed  into  a  varia¬ 
tional  principle  for  the  Lagrange  integral  functional  L  =  ///  L(r , t) d 3r ,  the 
extremals  of  which  give  the  Eulerian  equations  (1)  for  the  components  s  =  1,2,...N. 

In  this  approach,  Eqs.  (2)  -  (A)  are  auxiliary  equations  by  means  of  which  the 

^ 

fields  ng(r,t),  P  (r,t),  and  E(r,t)  can  be  eliminated  from  Eq.  (1).  This 
elimination  will  be  carried  through  after  introducing  canonical  variables. 

In  analogy  to  the  Lagrange  density  L(r,t)  for  nonconducting  fluids,—  the 
Lagrange  density  for  the  many-component  plasma  with  the  self-consistent  field 
E(r,t)  is  sought  in  the  form 


L( r,t) 


_>  2 
n  m  v  *■ 
s  s  s 


Y 


(7) 


If  Eq.  (7)  is  the  correct  Lagrange  density,  the  Eulerian  equations  of  motion 

6 

for  the  plasma  components  should  follow  from  the  extremum  condition  — 


<S/dt/ / / f (r,t)d3r  =  0, 


(8) 


r .  e. , 
t„ 


N 


fdtfff  )  [n  m  v  -dv  +  —  m  v?6n  - - - 

i  j  j  ‘  s  s  s  s  2  s  s  s  v  -  1 
t  2  s— 1  's 


6P 


-  e  E*  6E] d 3r  =  0. 
o 


(9) 


The  Eulerian  virtual  displacement,  Oy(r,t),  of  a  mass  element  of  the  s-th  com¬ 
ponent, 


► 

o 

s 


(r,L) 


v 

s 


(r, t)6t , 


(10) 


vanishes  at  the  initial  (t^)  and  final  ( 1 2 )  instants  of  the  .variation  of  the 
trajectory  of  the  fluid  element  which  occurs  in  the  interval  t^  <  t  <  t^» 


a  (r,t  )  =  5,  a  (r,t  )  =  5, 

s  1  s  £ 


(11) 
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and  is  tangential  to  tne  planes  dA  bounding  the  plasma  at  r  =  s : 


0s(r,t)  -dA(s)  =  0.  (12) 

The  variations  of  the  fields  of  each  component  s  are  given  as  functionals  of 
its  displacement  o^Cr.t): 


* 

iSv  - 

3  o  -»■■•>■  -»■ 

s  +  v  "S/a  -  a  *Vv  , 

.  ^  s  s  s  s 

(13) 

s 

3 1 

= 

-  V .  (n  o  )  , 

(14) 

s 

s  s 

6P  = 

-  o  -VP  -  y  P  V*o  , 

(13) 

s 

s  sss 

~y 

fN 

—  1  n  -V 

<5E  = 

-tinea 
o  ^ ,  sss 

(16) 

Equation  (13)  is  obtained  by  setting  the  variation  of  the  Lagrangian  velocity 
equal  to  the  substantial  derivative  of  the  Eulerian  displacement.  Equations 
(13)— (16)  result  directly  from  Eqs.  (2),  (3), and  (6),  by  means  of  Eq.  (10). 
Substitution  of  Eqs.  ( 1 4 ) — ( 1 5)  into  Eq.  (9)  gives 


-> 

Do 


/dt/jf  )  [n  m  v  •  (-r—  +  v  •  Vo  -  o  -Vv  )  -  ~r  m  vzV- (n  o  ) 
1  J  J  L ,  s  s  s  9 1  s  s  s  s  2  s  s  s  s 


t  U  s=  1 


+  — — : — (o  *VP  +  y  P  V*o  )  +  n  e  E-o  ]d3r  =  0 
Y  -  1  s  s  ss  s  ss  s 


(17) 


By  means  of  partial  integrations.  Gauss'  integral  theorem,  and  the  surface 
condition  (12),  Eq .  (17)  is  reduced  to 


t.. 


N 


fdtjjf  )  a  •  [—  (n  m  v  )  +  V*(n  mvv)+VP  -  ne  Eld3r  =  0  .  (18) 

.  K  3t  sss  ssss  s  ss 

t  SI  S=  1 


Since  Eq.  (18)  is  zero  for  arbitrary  variations  a  ,  s  =  1,2,...N,  the  integrand 
in  the  bracket  must  vanish  a‘:  any  point  (r,t)  of  the  plasma.  Thus,  the  Eulerian 


equations  of  motion  (l)  are  obtained  for  the  N  components  from  the  variational 


principle  (8).  This  completes  the  proof  that  liq .  (7)  is  the  correct  Lagrange 

density  for  the  multicomponent  plasma. 

The  Lagrange  density  L  and  the  Hamilton  density  If  are  not  only  functionals 

of  fields  of  r  and  t,  but  also  depend,  in  general,  explicitly  on  r  and  t.  In  the 

-y 

following  canonical  theory,  the  explicit  r,t  dependence  need  not  to  be  formally 
indicated,  since  the  virtual  displacements  constructing  the  varied  paths  are  for 


constant  r  and  t. 
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LAGRANGE  EQUATIONS 

!n  the  Lagrangian  formulation  of  many-component  plasma  dynamics,  the 

Lagrangian  L  =///  L(r,t)d3r  is  expressed  as  a  functional  of  the  fields  q^Cr.t), 

-y  ~y 

which  will  be  shown  to  form  with  the  fields  pg(r,t)  canonical  variables.  These 
are  defined  by,  s  =  1,2,...N: 

p  =  m  v  ,  3q  /3t  =  n  v  .  (19) 

s  s  s  s  s  s 

■+  ->  -y 

Note  that  pg  is  proportional  to  the  velocity  field  Vg,  whereas  qg  is  a  t-integral 
functional  of  the  particle  flux  nsvs>  Integration  of  the  auxilary  equations  (12), 
(3),  (6)  (Eq.  (6)  is  equivalent  to  Eq.  (4)],  and  Eq.  (19)  yield 


n  =  v  -  V  *q  ,  v  =  v  (r) 
s  s  s  s  s 


P  =  P  [(v  -  V* q  )/n  ] 

s  so  s  s  so 


>  ->  -1  -*• 

E  =  e  -  e  e  q  ,  e  =  e  (r) , 
S  OSS  s  s 


(20) 

(21) 

(22) 


v  =  (3q  / 3t)/(v  -  V-q  ) 

s  s  s  s 


(23) 


P  ,  n  ,  v  ,and  e  are  constants  with  respect  to  the  integration  variables. 

SO  so  s  s 

Substitution  of  Eqs.  (20)  -  (23)  into  Eq .  (7)  gives  the  Lagrange  density  of  the 


plasma  as  a  functional  of  the  qg, 
N  ^  ^  ^ 

L  =  I  LS(V  3qs/3t’  v'q  =  ) 


S=1 


(24) 


where 


i  =  %  m 
s  s 


( 3q  /3t)2  P 


v  -  V*q  3 

s  ^ s 


v  -  V  •  q  v  c  , 

so  ,  s  s  .  s  o  -I  t  v? 

- _  ( - ) - —  (E  _  c  e  q  )/ 

>  -  1  n  2soss 


(25) 


so 


is  the  Lagrange  density  of  the  s  component.  The  Lagrangians  of  the  plasma  and 


i 

'i 
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its  components  are  given  by 


•L  =  I  L  ,  L  =  J/J  id3?. 


(26) 


s=l 


The  fields  pg  and  defined  in  Eq.  (19)  are  canonical  variables  since  by 


Eqs.  (24)  -  (25) 
>  3L 


3q  /3t 
s 


S  3(3qg/3t)  vg  -  V-qg 


—  =  m  v 
s  s 


(27) 


Thus,  by  introducing  the  canonical  variable  qg,  JL  has  been  decomposed  into 

the  Lg  of  the  individual  components.  Note  that  the  original  Lagrange  density  in 

Eq .  (7)  does  not  permit  such  a  decomposition  due  to  the  coupling  of  the  components 

by  the  common,  nonlinear  field  energy  e  F.^/2. 

o 

For  the  Lagrange  density  of  the  form  of  F.q.  (24),  the  variational  principle 


in  Eq.  (8), 

r 


N  C2 


/  dt  Jff  Ld  3r  =  l  6  /  dt  ///  L  d3r  =  0, 


s=l  t„ 


(28) 


becomes,  since  the  Lg  are  independent  for  s  ~  l,2,...N, 
1  2 

<5  /  dt  [ff  L  (q  ,  3q  /3t,  V’q  )d3r  =  0. 

*  i s  s  s  s 

t  J  u 


Hence 

t.. 


3q  3  L 


f  dt  Jff  l~-5q  +  - - •6(3-t-%+  —V  6V-q  ]d3r  *  0 

t,  n  3q  3 (3q  /3t  3  V*  q 

1  s  s  s 


(29) 


(30) 


where 


6q  =  (3q  /3t)6t=nv6t=no  , 
s  s  s  s  s  s 


(31) 


<$q  (r,t.)  =  5,  5q  (r,t  )  =  5,  5q  (r,L)  -dA(s)  =  0, 
si  s  2  s  r=s 


(32) 
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by  Eqs.  (10)  -  (12).  Ry  means  of  the  partial  integrations, 

31  3q  31  t.,  C2  „  31 


l2  3 1  3q  ,  3 1  t  C2 

/ - -A-  -6  (^-)dt  =  [ - — - - 6q  1 

t  3(3qs/3t)  D(3qs/Dt)  t,  t 


-  /  3q_-  ± - dt, 

3(3qg/3t) 


1 


(33) 


ill 

S2  3V*q 


3L  >1 

— 3-  6  V  •  q  drr  =  i  - 2-  ;q 

D  V  ■  c] 


-  ///  aq0‘v( 


31 


Q 


s  3  V  •  q 


L-)d3r, 


(34) 


and  Eq .  (32),  the  extremum  principle  is  reduced  to 
1 2  3l  ,  31  31 

/dtj//l-^  -  £  -  v(~A.)].^  d3^  =  o 

t,  3q  3  ( a q  3t)  3V-q 

1  *s  S  S 


(35) 


The  condition  for  the  vanishing  of  this  integral  for  arbitrary  6qg  leads  to  the 
Lngrnngi an  equation  of  motion  of  the  s  component  of  the  plasma: 


3__ 
3  t 


31 


31 


31 


-T  +  V(— ^-)  =  , 


(36) 


3(3qg/Dt)  3qg  37-q^ 


whore 


3q  /3t 


w  =  m  (- 
s  s 


)  X  V  x  (- 


3q  /3t 


(37) 


V  -  V'q 

s  s 


v  -  V*  q 
s  's 


is  a  vector  field  perpendicular  to  6q  ,  w  ■  5q  =  0  by  Eq.  (31).  It  is  noted  that 
the  extremals  of  Eq .  (35)  are  undetermined  up  to  fields 

>  - >  -*  >  • r  .  -  >  ■>  . 

w  i  iSii  (w  .  m  v  x  V  y  v  7  0  for  v  f  -Vi p  ). 

s  ssss  s  s  s 


In  terms  of  the  I.agrangian  of  the  multi-component  plasma,  Eq .  (26),  and 
9 

Junctional  derivatives,—  the  Lagrange  equation  (36)  assumes  the  most  general 
form, 

3  jI,  61. 


1  6(3q  /3t  5q 

s  £ 


(38) 


where 


41, 

-> 

5q.. 


31 

-r-7(— -J- 

3,,a 


6L 


31 


3  V  •  q 


6(3q  /3t)  3(3 q  /3t) 

S  S 


(39) 
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Subst  i  t  ut  ion  of  L  f  tom  Fq.  (25)  into  Kq .  (36)  indeed  yields  the  fc 


iquation  (l). 


.s  into  V  ( v  1 ) 
s 


^  -  + 

v  *  Vv 
s  s 


i 


1 e  r  Ian 
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HAMILTON  EQUATIONS 


Another  important  formulation  ot  the  dynamics  of  a  many-component  plasma  is 


obtained  by  means  of  the  Hamilton  forma  1  i  ;;m. 
momenta,  the  Hamilton  density  of  a  component  s  is 
3q 

H  =  p  ~  -  L  , 

s  s  Dt  s 


Since  p  and  q  are  canonical 
s  s 


.  p  v  V  •  q  ,  c  ,  „ 

H  =  - —  (v  -  v-q  )  P^  +  - 7  ( - )  +  -7  (t  -  e  e  q  )  (41) 

s  2m  s  s  s  Y  -  1  n  2  s  oss 

s  s  so 

by  Eq.  (25).  Accordingly,  the  Hamilton  density  of  the  plasma  is  a  functional  of 
the  form 

N 

H  =  'i  H  (p  ,  q  ,  V-q  )  .  (42) 

u  s  s  s  s 

S=1 

The  resulting  Hamiltonians  of  the  plasma  and.  Its  components  are 

H  =  [  H  H  =  ///  H  d3r.  (41) 

s— 1  b  s  n 

Thus,  also  within  the  Hamilton  approach,  a  decomposition  of  the  Hamiltonian 

H  into  component  Hamiltonians  H  is  achieved,  each  11  depending  only  on  the 

s  s 

fields  of  the  component  s.  This  is  quite  remarkable,  since  the  various  components 

r  > 

s  Interact  through  the  self-consistent  field  E(r,t).  The  decomposition  is 
physically  possible,  since  the  electric  fields  of  the  different  components 
superimpose  in  a  linear  way. 

The  Hamilton  equations  of  motion  of  the  components  s  follow  from  the  action 
principle  (29), 


<5  /  dt  ///  (P/T,:*  -  =  °. 


(44) 


which  gives,  since  the  Hamilton  density  is  a  function  of  the  form 


Hs =  Hs(lv  v  V‘i)’ 

*~2  Dq  3H  ^  3q  a  1/  jH 

!  dt  IS!  Ujt- — r>'6Pa  +  Ps'6  (  3T  ) - - i“  * '■’*%] d  r  =  0  .  (45) 

t,  ft  3p  ‘  oq  jV • q 

I  s  s  ns 

Interchange  of  the  sequence  of  variation  and  differentiation  and  partial  inte¬ 


grations  , 


2,  -  >  3 


t., 

t0  l  .v- 

2  f  ox 


/?8*FqsdtBiv\it  -  ri±  -*vit. 


li  ■  1  4  :)l 

3tf  "jH  dH 

US  - -  V-<5q\d3r  =  if - ~  <Sq  *dA  -  ///  6q  -V( - -)d3r,  (47) 

ft  3V.q.  b  3V-q,  S  ft  S  3V-q 

b  S  S 

where  t he  variation  5q^  satisfies  the  condition  in  Eq,  (32),  transform  Eq.  (45) 


t2  9q  3  H  3p  3  H  3  H 

1  dtj//  ~  ~)’6p  -  +  V(— S-)].«q  >d*r  =  0.  (48; 

C1  «  3PS  3^s  3^s 

The  conditions  for  the  vanishing  of  this  integral  for  arbitrary  variations 

->  -V 

4p^  and  5q^  gives  the  Hamilton  equations  of  motion  for  the  plasma  components  s: 

3q  3H 

s  s 

IT =  ~  .  (Ay; 

3ps 

3p  3 II  '  3  H 

s  s  „ ,  s  .  -> 

*T  +—  -  v(— — )  =  wg  ,  (50) 

3q  __  d7*q 


where 


->  1  >■  -> 
w  =  —  p  x  V  x  p 

*■’  TTI  ft  ft 


m  s 
s 


is  a  vector  field  perpendicular  to  iq  =  n  v  o t ,  i . e . , w  • 6q  =0.  As  above 

s  s  s  s  s  ’ 

ws  is  determined  from  the  condition  of  rotational  flow  (w  f  ft  for  V  x  v  /  S). 
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By  introducing  the  Hamiltonian  ni  the  p ! asina  ,  Eq .  (ah),  and  functional 

Q 

derivatives,-  the  Hamilton  equations  (49)  -  (50)  are  brought  into  the  most 
general  form. 


at 


> 


(52) 


3t 


4H 


(53) 


where 


6H 

* 

OH 

s 

-►  ’ 

V  jHs 

Oil  b 

51/ 

-  V  (— J) 

(54) 

"Ps 

Dp 

s 

fiq  Dq 

s  s 

3V-qs 

Evaluation  of  Eqs.  (49)  -  (50)  for  the  Hamilton  density  in  Eq.  (41)  results  in 
Eq.  (19)  and  the  Eulerian  equation  of  motion  (1),  respectively. 


CONCLUSIONS 


Comprehensive  Lagrangian  and  Hamiltonian  methods  have  been  developed  for 

many-component  plasmas  with  internal  interactions  through  the  self-consistent 

electric  field.  The  theory  is  based  on  the  Lagrange  density  (26)  and  the 

Hamilton  density  (41),  and  the  canonical  variables  p^  and  q  defined  by 
'  ►  ->  .  '► 

p  =  in  v  and  3q  / D t  =  n  v  .  respectively.  In  the  limit  of  a  single  component 
b  s  s  s  s  s 

(N=l)  and  vanishing  Coulomb  interactions,  the  corresponding  Lagrange  and 

Hamilton  equations  of  motion  for  an  ordinary  neutral  particle  gas  result. 

The  canonical  theory  presented  reduces  the  coupled  partial  differential 

equations  (1)  -  (4)  for  the  fields  v  (r,t),  n  (r,t),  P  (r,t),  and  E(r,t)  to 

s  s  s 

a  single  partial  differential  equation  for  the  canonical  field  qg(r,t).  The 
fundamental  Lagrange  equation  of  motion  is 
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mathematical  advantages  of  the  theory  presented  are  obvious,  in  particular  for 

plasmas  with  several  components.  In  comparison  with  the  previous  formulations 

13  4  5 

of  variational  principles  in  fluid  mechanics  — * — * — *—  which  make  use  of  an 


excessive  number  of  scalar  Clebsch  parameters,  the  present  canonical  theory  in 

-+  -> 

terms  of  the  canonical  vector  fields  pg  and  excells  in  mathematical  simplicity. 
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VI.  Distribution  function  of  turbulent  velocities 


ABSTRACT 

The  Hamiltonian  H  and  Hamilton  equations  of  motion  are  derived  for 
-►  -► 

the  Fourier  amplitudes  p£  and  q+  of  the  canonical  conjugate  fields  p  and 
•+  ->  ■+  “►  -► 

q  defined  p  =  mv  and  Dq/3t  *  nv,  vhere  v(r,t)  and  n(r,t)  are  the  velocity 
and  density  fields  of  ideal,  compressible  gases  in  the  state  of  fully 
developed  turbulence.  A  Liouville  equation  is  presented  for  the 
distribution  function  f (p^,q^; {lc})  in  the  multidimensional  phase  space 
formed  by  the  scalar  components  of  the  set  {£}  of  wave  mode  vectors 
P£  and  q^.  As  an  application,  that  stationary  solution  f  «  f (It)  of  the 
Liouville  equation  is  calculated,  which  maximizes  the  turbulence  entropy. 

It  is  shown  that  the  distribution  of  the  velocity  and  density  fluctuations 
of  compressible  gases  is  Gaussian  in  fully  developed  turbulence,  in 
agreement  with  the  experiments. 


INTRODUCTION 


The  turbulence  probLem  has  aroused  the  interest  of  many  celebrated 

1)  2)  i)  4) 

researchers,  such  as  Reynolds  ,  von  K Arman  ,  Taylor' ',  Heisenberg  , 

Kolmogorov  ,  Prandtl^ ,  and  Frenkiel ^ ,  to  name  only  a  few.  In  spite 

r  ,,  ,  8-10)  ,  L 

of  all  progress  made  to  date  ,  many  experimental  observations  on  the 
turbulence  phenomenon  are  still  unexplained  from  first  principles.  One 

of  the  most  elementary  experimental  results  is  the  Gaussian  velocity 

-►2  -2  -> 
distribution  ~  exp(-v  /c  )  of  the  turbulent  velocity  fluctuation  v  in 

11 )  * 

compressible  gases.  The  Gaussian  velocity  distribution  of  fully  developed 

turbulence  is  a  fundamental  fact  as  the  Maxwellian  distribution 

1  ->2  ■> 

~  exp(— ^mv  /KT)  of  the  molecular  velocities  v  in  statistical  equilibrium 

gases . 

The  Maxwellian  velocity  distribution,  which  is  the  result  of  many 
molecular  collisions,  is  independent  of  the  molecular  interactions,  which 
have  brought  about  the  statistical  equilibrium.  The  experiments  indicate 
that  in  fully  developed  turbulence,  the  distribution  function  of  the 
velocity  fluctuations  shows  only  in  the  far-out  wings  jv  ]-*•  «  smal  1  deviations 
from  the  Gaussian  distribution . ^ ^  Apparently,  the  distribution  of  the 
velocity  fluctions  in  stationary  turbulence  depends  hardly  on  the  viscous 
momentum  transfer,  which  is  negligible  for  the  velocity  fluctuations  of 
large  scale  but  significant  for  the  velocity  fluctuations  of  small 

scale  A->  0  ,  concerning  the  sustainment  of  the  turbulent  state  since 

-V^v,*-  ~  (A  =  2-n/k)  . 

k  k 

In  other  words,  the  Gaussian  distribution  of  the  velocity  fluctuations 
in  fully  developed  turbuLence  is  mainly  sustained  by  nonlinear  mode  coupling 
of  the  turbulence  elements  of  different  wave  vectors  ic.  For  this  reason. 


it  should  be  possible  to  derive  the  turbulent  velocity  distribution  from 
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a  statistical  hydrodynamics  of  the  turbulent  llnid  notions,  without  taking 
viscosity  into  account.  By  Fourier  analyzing  the  canonical  conjugate 

-  V  -  ►  -+ 

variables  p(r,t)  and  q(r,t)  of  the  turbulent  gas,  it  is  shown  that  the 
->  -> 

Fourier  amplitudes  p^(t)  and  q£(t)  satisfy  Hamilton  equations  of  motion,  i.e. 
the  random  wave  modes  k  behave  like  quasi-particles.  This  makes  it  possible 
to  derive  a  Liouville  equation  in  the  hyper-space  lormed  by  the  scalar 

•  V  -r 

components  of  the  set  lp^  ,  q^}  .  From  this  fundamental  equation,  the 
quasi -Gaussian  distribution  for  the  turbulent  velocity  fluctuations  is 
derived  as  that  stationary  solution  which  maximizes  l  he-  entropy  of 


turbulence.! . 


HAMILTON  EQUATIONS 


In  the  statistical  turbulence  theory  to  he  presented,  tiie  effects  of 
viscosity  are  disregarded  so  that  the  random  velocity  field  v(r,t)  of  the  gas 

-  t-  —V 

Is  i rrotational ,  Vxv=0,  in  the  absence  of  external  forces  which  could 

>  > 

artifically  produce  velocity  curls.  The  condition  Vxv=0  defines  longitudinal 

or  so-called  "acoustic"  turbulence  in  compressible  gases. It  represents 

a  good  approximation  for  the  treatment  of  certain  turbulence  phenomenon, 

such  as  the  present  calculate  on  of  the  turbulent  velocity  distribution. 

The  basis  for  the  continuum-statistical  considerations  form  the  Hamilton 

equations  of  motion  for  the  compressible,  nonviscous  gas.  If  v(r,t) 

•  ► 

and  n(r,t)  designate  the  velocity  and  density  field  ol  the  gas,  respectively, 

■v  >  ir  ■  v 

the  canonical  conjugate  fields  p(r,t)  and  q(r,t)  of  the  gas  are  defined  by  the 
proportionality  (m=mass  of  gas  molecule)  and  functional  relations , respectively 


p  =  mv 


3q/3t  =  nv 


The  Hamilton  equations  of  motion  and  the  Hamiltonian  H  of  the  gaseous  system 

of  volume  ft  =  f f f  d  *r  =  L  L  L  are  given  by : * ^ 

x  y  z 

3q  _  Mi 

3jV  =  _  6H 

3t  6cf  ,  (3) 


,  1  V  ^  ()  /  v~  V  * 

H  =  -M’q)p  +  — f  — 


=  3H 

6  p  3  p 


fib  =  _3H  _  „  /_QH_ > 
6  3  ~  3d  '  \3V-q 


delines  the  functional  detivatives  4  |3H/3q  0  since  If  =  H(p,V-q)  by  Eq.(A)|. 


In  Kq.(A),  y  (,  =  (■/<•)  is  Lite  aiMabati,-  coefficient,  P  and  n 

o  o  1  V  o  o 
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[reference  values  of  P ( r , L )  and  n(r,t)]  and  v  are  constants  obtained  through 
integration  of  the  polytropic  energy  conservation  and  continuity  equations  of  the 
gas,  respectively. 

Since  n(r,t)  and  v(r,t)  are  random  turbulent  fields,  p(r,t)  and  q(r,t) 
are  represented  as  complex  Fourier  series, 

P  ( r ,  t )  =£Pjj(t)  exp  (ik-r),  q(r,t)  =^q^(t)  exp  (ik-r),  (6) 


where 

k  =  2tt 


V  M  \>  I 

x  _ z\ 


L  ’  1.  L 
x  y  z 


{ k} 

=  0,  ±  1,  +  2, 


,  ±  v 


(7) 


x,y  ,2 

are  the  discrete  eigenvalues  (Q<°>)  of  the  wave  vector  k.  Since  the  minimum 
wave  length  X  =  2n/k^  is  of  the  order  of  the  mean  free  path,  Z,  one  has 
v  ~[L/£]>>  1  for  Z<< L. 

Hamiltonian  equations  of  motion  similar  to  those  for  the  canonical- 

■  V  V  -k  -  F 

conjugate  fields  p(r,t)  and  q(r,t)  can  be  derived  for  their  Fourier  amplitudes 

V  — V 

p,>(t)  and  q,k(t),  i.e.,  for  the  individual  wave  modes.  By  Eqs .  (2)-(b), 

q,i/ dt  (8) 


311 


Dji 


-k 


t)  ,  i.e.,  tor  the  individual  wave  modes.  By  Lqs . 

for 


!l 


52 


and 


4  ■  f  te 


57 


3  H 


3? 


+v.  Mi.  = 


I- 


ik'rV'[-Jl  d3 


-  52dp//dt 


(9) 


<  ►  3  -v  ij< .  r 

since  Jj/7- A  d  r  =  0  if  A  is  expandable  in  orthogonal  functions  c  =  unit 


tensor) . 

The  canonical  equations  for  the  canonical -conjugate  amplitudes  p^v C t )  and 
qj>(t)  of  the  individual  wave  modes  k  are  by  Eqs.  (8)-(9): 
dP2 


1  3H 


dt 


51  3<t 


(10) 


-k 


and 


dq. 


1  311 


dt  J 


(ID 


-k 


75 


Equations  (10)  and  (11)  represent  a  complete  system  of  coupled  nonlinear 
differentia]  equations,  which  are  antisymmetric  in  the  k  indices.  These 
canonical  equations  determine  the  temporal  development  p'(t)  and  q  +  (t)  of 

K  K. 

>■  -  r 

-►  i  k  *  r 

the  wave  modes  k  of  spatial  structure  =  e 


In  Eqs.  (10)— (11)  the  canonical  equations  of  the  wave  mode  dynamics  are 
stated  in  complex  form  for  the  conjugate  variables 


k  = 

k*  k =  real 

(12) 

->  -V 

-y 

Si 

1! 

Qk  +  1S'£  * 

Qk%  S-  =  real  , 

(13) 

so  that 

> 

-y  .  *v  -y 

Pk  = 

2(pk  +  p-k} * 

k  =  0/2i)  (p£_p_£)  » 

(14) 

2(qk  +  q-k) * 

-y  ~y  -  *■ 

^k  = 

S£  =  (l/2i)  (qk-q_k)  , 

(15) 

since 

-y 

—rh  -y 

(16) 

P£  = 

p-k’ 

qk  = 

for  real 

vector  fields 

p(r,t)  and  q(r,t),  Eq.(6).  Substitution 

of  Eqs . (12) 

(13)  into 

Eqs.  (10)-(11)  yields  the  Hamilton  equations  for  the 

real  and 

imaginary 

parts  of  the 

canonical  conjugate  mode  amplitudes: 

dP,>- 

k 

1  3H 

dPk  1  Dll 

(17) 

dt 

"  '  n  ’ 

dt  Q  3?- 

k 

d% 

1  3H 

d?5>-  ,  ... 

dt 

k 

k  _  1  311 

dt  SI  3l+ 

(18) 

li 


MOUVILEE  EQUATION 


The  existence  of  the  canonical  equations  (10)  and  (11),  which  arc- 
analogous  to  the  canonical  equations  of  a  many-hody  system,  indicates  a 
quasiparticle  behavior  of  the  individual  wave  modes.  For  statistical  purposes, 
let  a  large  number  (ensemble)  of  similar,  turbulent  gaseous  systems  be 
introduced  (each  containing  a  set  { k  1  of  wave  modes  or  quas  ipartic  1  es  Ic) 
in  the  mul  tidemensional  phase  space  formed  l.y  the  real,  and  imaginary  parts 
of  the  Fourier  vector  components  p^>  .  and  q^  ^  (i  =  l,2,  1)  of  the  entire  set 
ikjdefined  in  Eq.(7).  The  number  of  gas  systems,  which  have  their  phases  in 
the  volume  element  il  ^  ^  (d  p^d^q  j>)  at  the  point  (p£,  q^;{k})  of  the  phase 
space,  is  given  in  terms  of  the  density  in  phase  space,  f=f(p>,  qj>,  t;  {It}),  by 

K  K 


dw  =  t'(p£,  qkr,  t;  >  k  I  )  H  (d  pod  q^)  , 

{ k } 


whe  re 


‘  dV’Ri ;•  ■  dV%-  (20) 

in  accordance  with  Eqs .  (12)  and  (13).  The  phase-space  density  satisfies  the 
continuity  equation  in  phase  space,  since  the  number  cj f  plasma  systems  in  the 
ensemble  is  conserved. 


Ft  +  "(it)  '  <v(k>£)  *  0  • 


where 


.vfei  ^ 

■  kl  Z-  l<it  '  dt 


\ 

«> =  IrP* ' 


The  canonical  equations  (LO)  and  (11),  which  determine  the  motion  of  any 


point  in  phase  space,  indicate  that  the  phase-space  fluid  is  incompressible. 
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;►  1  X — ' 

7{k}‘  ik}  ■  ri  ]C 

(ki 

because  ol  this  property,  the  continuity  < •  ( | u 1 1  inn  (21)  can  he  brought  into  the 
form  of  the  so-called  Liouvi-llo  equation: 


This  fundamental  equation  describes  the  temporal  development  of  the  distribution 
in  phase  space.  ft  extends  statistical  mechanics  to  gas  continua  .  in 
random  motion. 

The  Hamiltonian  U  in  Eq .  (4)  is  an  invariant  for  a  nondissipative  gas. 

The  steady-state  solutions  of  Kq .  (2b)  are  functions  of  H : 

f  =  f  00  ,  (26) 


since 


21 

dt 


I  V'/lL  >->H _ ah  iLL  \  af 


0,  for  f  =  f (II) 


(27) 


In  applications,  the  compact  complex  firm  (2b)  of  the  hiouville  equation 
is  mathematically  preferable.  In  terms  ot  real  vector  components,  Kqs.  (22)- 
(23)  become 

bit  ■  E 

{ k> 


f^'j 

l  d  t 


dR,* 

_  k 

dt 


d^k 

dt 


(IS; 


dt 


(28) 


(  k) 


(29) 


Since  the  phase  space  fluid  is  i ncompressib I e , V ^  j  * 

suit  st  i  t  lit  ion  of  Kqs  .  (28)- (29 )  into  the  continuity  equation 

the  hiouville  equation  in  real  notation: 

a  i_  .» n _  b[  _  ri,i_ 

Vnk  '  T*k  '  ^  ^ 


=  0  by  Eq.  (24) , 
(21)  results  in 


1 


7K 


Tl ic*  solution  of  the  Louville  equation  permits  the  determination  of  tli 
macroscopic  properties  of  turbulent  gases  as  ensemble  averages.  The 
statistical  considerations  are  based  on  the  quas i -part ic ] e  character  of  th 
individual  wave  modes. 
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TllRBlII.KNCK  MNTROPY 


The  state  (if  stationary,  fully  developed  turbulence  is  one  of  maximum 
probability  W.  Since  the  entropy  of  a  turbulent  system  is  a  function 
S  "  fn  W  of  its  probability  U',  the  entropy  oi  turbulence  must,  assume 
an  extremum  in  stationary  state. 

The  entropy  associated  with  the  turbulent  modes  of  the  gas  is  defined 

-y  -y  y 

as  the  phase  space  integral  over  the  distribution  function  f  =  f (p^,q-+; (k)  ) 

K  K. 


■H 


S  =  -K  |  ...  |f  fllf  H  (d6p->dfVq>) 

{£}  k  k 


('ll) 


where  k  is  a  dimensional  constant.  This  definition  of  turbulence  entropy  is 
in  complete  analogy  of  the  definition  of  turbulence  for  arbitrary  nonoqui li¬ 
brium  systems. 

The  idealized  description  of  the  gas  as  a  continuum  without  collisional 
dissipation  is  reflected  in  the  time  -  indepencence  of  the  turbulence  entropy 
By  I’qs.  (3  ),  (21),  and  (24). 


dS 

dt 


cj  +  Inf)  ^  +  v 


<3t 


(kl 


v,-:-,  ;f  n  ( 

lk}J  Ik) 


,6- 


f  ydVdV)  =  0.  (32). 
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VELOCITY  DISTRIBUTION 

As  an  application  of  the  canonical  theory,  the  statistical  distribution 

of  the  random  velocities  in  turbulent,  compressible  gases  is  calculated. 

Let  the  turbulent  pas  be,  on  Lhe  average,  homogeneous  (n  =■  n  )  and  at  rest 

o 

(v  -  0)  so  that  for  the  spatial  mean  values 

p  =  5  ,  q  +  d  .  03) 

The  gas  is  assumed  to  be  in  contact  with  an  appropriate  other  system  which 
sustains  stationary,  fully  developed  turbulence.  The  canonical  conjugate 


fluctuation  fields  can  then  be  expanded  in  the  stationary  Fourier  series: 

p(t)  =  2_J  Pk*f  txplik-r  I,  q(r)  =  2_j  q^exp  [  ik  •  r  ] ,  (14) 

{k}  f  k  > 

where 


10) 


Pk'  =  ^  j"j'j*p<r)ox|»[-iE-r]d  (r  , 


q(r)exp[-ik-r]d  r 


(33) 

(3b) 


are  the  complex,  time-independent  Fourier  amplitudes  of  the  real  fields  p(r) 
and  q (r) : 


•>*  _  -* 

PC  =  P-k’ 


-V 

qk  =  q-k 


CD 


The  orthogonality  relations  for  elementary  vector  fields  of  the  form  given  in 
the  above  Fourier  expansions  are 

JJJ(P  )k‘  (q)k,d  r  =  !JJPk 'qlt'oxpfl(k  “  k'),rld'r  =  u  -CD 

Hy  Eq .(27),  the  phase  space  density  of  a  conservative  gas  (U=E)  in  a 
stationary  state  in  on  arbitrary  function  of  the  Hamiltonian 


.  q£,  f k 1  )  =  f[II(p£,  q+;  Ik} 


09) 


a  i  ua  t  state,  tiie  entropy  of  turbulence  assumes  an  extremum (state 
r  .  b v  l.q  .  (31  ) 


HI 


S(P£,qk-;  ik  I 


)  =  /r<  Pk‘,qs  fk)  )  e»  i(Pk-.^;  no  )  :i(d  V  V‘k 

1  k  i 


)  .  (40) 


Ai  cord  inn !  y  >  llie  phase-space  d  ist  r  iliut  ion  1(11)  ol  .stationary,  lullv  deve  1  oped 
turhuU’ine,  is  determined  by  the  variational  equation 

J.  "*  <’ti  O  li  (d  |),-d'  <!,• I  0  ,  (41) 

where 


-  -  k  I.  .  .  I  (  1  I  t’li  1)  sf  11  (d  ^|i  -df  (|  •  I  ■> 

i  k  } 


J...Jilll(d,,1>'t;d,,q^)  =  0  , 

f.  .  .  f  H  M  i:  (d'V/^q,  )  -  0  , 

3  3  \i\  k  " 


t  > 


f  S) 


since  tin1  number  of  jrik  systems  and  the  tot  il  <  ■  1 1  •  rp.y  ire  i .  <  1 1  i .  u .  *  .  •  ‘  |  - 

ii.u  1 1  i  p  I  i  * -.it  i  mi  i  •  1  K*| .  ( 4 1’ )  by  (  <-  1  1  and  Uq  .  <  4  ’■ )  l>  >  !■  ,  ■  '  <  ■  i  •  !  in  <  •  l 
iimJi  t  *■  mi i  nod  kayraa;',  i  an  mu  1 1  ip  I  i>  r.-,  i  v* •  ■ 

<■  r  I  ■  i-.  • 

J.  .  .J  <  t  W  t  +  I  +  .  It)  l  _U!  p^-d  1^1  II 


1  +  *  ■*  I!  -  0 


i  s 

1  1  '  '  i  1  I  ■  !  1  V 

.  l-'i  'in  !  <|  . 

4  ■  )  tin 

rial 

i  .’.it  i  i'ii 

(’ll 

[-.H|-  (d*! 

•  k  :• 

t>  - 

k’1’  V 

in  t  1  ■  t  i u  ni 


f  ( H  > 


■  k  ■■ )  I 

/f  7  .  ,  4  •  fi  > 

•  •  -J  I  -  • 4i(fv,<|;;  Ik.  <)  }  [  (d  p-d  q^-). 

k  1 


( 4 ;  i 


[he  Hami  1  toil  i  an  in  i  '.<  j  .  (4  )  becomes  alt*  r  substitut  ion  **1  'be  expansion 

i  n  Kqs  .  ()’>)  —  (  !b  ) 

II 


1 4  H ) 


-  ‘  ^k-q  >exp|  i 


n  "TT,  k 
< )  \  K  I 


-  k  •  1 1 


.1  r 


( 49 ) 


-M 


since  v 
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=  n^.  An  evident  bionomial  expansion  and  subsequent  spatial  integration 
gives,  under  consideration  of  Eq .  (38), 


where 


II  -  i: [-  — ar —  +  ~  E  pep,*  +  1-yl’  n  2  E  (ik'd 

L'"1  2m  un  k  k  2  °°  :Ki 


k)(,k-q^) 


f  h(pk,qk’<k  ^50) 


.-1 


c-o' 


-J  W  n  .'(1)  .*(/-]), 

t  U  l 


E  (ikn)-qy ° i...ikv  ’ f i , 1  j; it 


:('-i )■'  (/-]) 

'k 


|kU)  +  ...k(  "1)l-q_(k-(l)  +  _  _  .  .-(''-!) 


! 


■(kv"  +  ...  +  k'  ij 
contains  the  terms  ot  third  and  higher  order  in  tin  Fourier  amplitudes.  The 

various  sums  in  Eq .  (3D)  have  the  physical  •  -•  *. ,  i»  i  np.  el  spatial  averages 

( t1.  a  e  v.  i  I  1  s  t  lien  I  a  in) 


151) 


V 


Vk  =  P1  ,:|  v 


'52) 


(  i k  ’<)/)  (  i  i  "l|”  1  E'V  "k*  5 

ik  1 


i  >  1 1 


win  i  » ■  v  (v  I')  and  u  -  n  -  u  ate  the  :  >  h  ■  1 1  <  •  t  •  t  veioc  it','  util  tlete.il'. 

t  I  it,  I  i  la  t  i  or.  s  a  t  the  comp  less  i  h  1  e  re.. 

<‘e"ih  i  tiitij'  el  la|s.  (\i)  and  (Mi)  vield-  the  .t  it  inuarv  i  .  t  r  i  b  u  t  i  on  1  ii:i<  t  im: 

in  tin-  to  mu  : 


t  (I!  ) 


■  to  mu  :  r  _ _ „ 

I  ;’oV  •  >  *  1  -  ’ V  •  • 

•  U..7fk'Pk  f  2 '  :’ono  ^Mk-qk 


I  (  i  k  -  q,*  ) 


-.<h<Pyqk*:<k*})J 


.  >'*  ) 


I  e  >••!>(;•■  I'  (<-I)  1  I  J.  .  expl-l  ll(p^  ,q,-;  k  )]  il  dhp^df*q^ 

i  k  i 


!,/  introducing  the  spatial  averages  defined  in  laps.  (52)-(3  3)  a  physically 
more  i  I  1  ust  rat  i  ve  iopresfnt.it  inn  <>f  the  stationary  turbulence  distribution 
i  s  o!i  t  a  i  tied  : 


I  (II) 


A  expj-  r'..jj,n(|m  v2  +  yyl’ 0(n/n(1) 2J  -  Fh(p^.q^;  (k  t) 


f 


(55) 


wile’ re 


]  =  [oxpU-^JMy -I)  1 )  1 J*. .  .Jexp[-3H(p  *,qg;  (k  })  )  I!  d6podf'q^ 


The  parameter  3  of  the  distribution,  which  has  the  dimension  of  a  reciprocal 
energy,  is  given  in  terms  of  the  total  r  i  i  i  l>n  1  <  n  t  energy  K  contained  in 
the  I  I  net unt inns  through  Lhe  relation 

p 

r  /...j*  IKp^c^tlf ')!  (plKpgoJ-tlkl)  U  d<>|»|_-ii,’«I‘k-  -  ■  (  ,  (  ,1 ) 

Ik1 

in  v,  i  ew  c>f  the  noi  m  a  1  i  xat  i  on  delined  in  Iq.  !\h).  f  /  <  ,  —  I  j  represents  the 

o 

non  t  u  rhu  1  eii  L  or  thet.;  il  energy  of  tin-  gaseous  system  of  >.  o  1  uii’e 

la|.(Si)  indicates  that  t  he  d  i  s  l  r  ilnit  i  on  I  n.i.  t  ,o,i  f  ( M  )  lor  the  turhulent 
1  1  in  t  u.it  ions  lias  the  tallowing  prnpertii 


it  Tlie  d  i  s  t  r  1 1  ut  i  on  .  >  1  t  he  v >  ■  I  o,  i  t  ,  ‘  I  s.  t  u.  ■  t  i  oie  i 


"  i  an  , 


.  ,  nir  •  (  jH) 

I  ( v  )  c  1  o 

i  i  )  Tin’  d  i  s  t  r  i  bu  t  i  on  o !  the  density  :  I  u.  t  ua  t  i  oils  is  o.  cu..  .  i  an, 

.  , '  2.  -\,V  n  .  .«*’  .  (V)> 

t  ( n  >  c-  2  o  o 

ihe  h  i  ghc  r  ctchr  t  c  mis  in  the  I  i  s  t  i  i  l.ut  i  mi  I  nm  t  I  on  i  f  1 1  )  -  •  I  I'q .  ( Vi  ;  I  ue 


to  h  (  p  •  ,c|  • ,  ;  k  • )  ,  la  | .  ( '»  I  )  .  a  t  e  ditti  nit  l 

e  K 

lot  low  i  tig  reasons  .  I  In  v  represent  <  i 


i  vi'  in  i  :  i  r  intent  •.  I  m  tli 


•  ■  I  s  o  1  eh  -ns  ;  t  s  aid  a-  1  :s ■  i  t 


t  1ii>  tii  it  ic  is,  i  . .  .  they  are  uiiohsc  i  vab  ’  •  •  it  tin  ,  /pet  '  nr  a  1  he  ■  a  '  >r  r  •  its! 
only  velocity  or  dens  i  l  v  I  !m  t  it  it  ions.  'In.  i  ,  i  m  h  ( t  ^  ,  c| ,  - ;  .  I  )  i . .  snail  el 
third  and  higher  odor  in  the  f  1  in  t  u.it  i  or  amplitude;  and  is,  in  general, 
significantly  smaller  than  the  leading  second  order  terms  in  the  expansion 
"t  H<pk>,>ik-; 'k  1)  . 

Normal  distributions  of  the  forms  given  in  1‘cjs.  (riH )-('>'))  have  been 
observed  in  fully  developed  l  m  i>;  1  ence  of  compressible  gases.' Kr.per  imental  1  y 
the  distribution  ol  turbulent  velocities  is  obtained  by  plotting,  the  I requency 
of  the  occurence  (over  a  larger  time  period)  of  a  velocity  signal  with 
amplitudes  between  |v|  and  |v|  +  j  Av |  vs.  |v | ,  i.e.  by  a  time-averaging  process. 


Ill  the  distribution  function  of  Eq .  (55),  the  quantities  v 


and  n 

represent  spatial  averages  (Eqs.  (52)  and  (51).  This  formal  difference 
is,  of  course,  no  contradiction  between  experiment  and  theory,  since  the 
time  average  is  equal  to  the  spatial  (ensemble)  average  in  fully 
developed,  homogeneous  turbulence.  In  Fig.  1,  the  experimental  (o)  and  theo¬ 
retical  (-)  distributions  of  the  turbulent  velocities  ire  compared!'^ 

It  is  seen  that  the  theoretical  (emission  distribution  of  the  turbulent 
velocities  !s  in  excellent  agreement  with  tie  e.-.per  imental  dale. 

The  theory  presented  is  highly  idealized  since  viscous  and  thermal 
dissipation  in  t  iie  turbulent  gas  .ire  not  taken,  explicitly,  into  account  . 
it  represents,  however,  a  first  attempt  at  extending  statistical  mechanics 
to  r.indori  continuous  media,  such  as  turbulent  gases.  Extensions  of  the 
theory  to  include  viscous  and  thermal  dissipation  represent  mainly 
mat  hemal  i cal  problems. 


i 
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VII.  Appendix:  stress  relaxation  naves  in  fluids 


ABSTRACT 

The  Navler-Stokes  equations  for  incompressible  and  compressible  fluids 
are  generalized  by  inclusion  of  viscous  stress  relaxation,  as  required  by 
kinetic  theory.  Two  initial-boundary-value  problems  of  the  nonlinear,  generalized 
Nav ier-Stokes  equations  are  solved  analytically,  which  describe  the  propagation 
of  transverse  or  shear  waves  due  to  temporal  and  spatial  velocity  pulses  v(0,t) 
and  v(x,0),  respectively.  It  is  shown  that  transverse  perturbations  propagate 
in  form  of  a  discontinuous  wave  with  a  finite  wave  speed  due  to  viscous  stress 
relaxation,  whereas  the  conventional  Navler-Stokes  equations  result  in  non¬ 
physical  solutions  suggesting  a  diffusion  process  covering  the  entire  fluid 
with  ml i nile  speed. 


! 


INTRODUCTION 


The  nonlinear  incompressible  and  compressible  Navier-Slokes  equations 
represent  (quasi)  parabolic  and  hyperbolic,  partial  differential  equations, 
respectively.  The  former  propagate  signals  with  infinite  speed  and  the  latter 
propagate  certain  signals  with  finite  speed  in  fluids.  In  an  infinite,  homo¬ 
geneous  fluid,  consider  a  small  (linear)  velocity  perturbation,  which  is  re¬ 
presentable  as  the  Fourier  integral, 

-f-  ->  -* 

'*/  f  ►/IN  iwt  -  ik*r  ,  t 

v(r,t)  -  J  v(k)  e  aw  * 

—  oo 

over  elementary  waves  of  wave  length  >  =  2 tt / k  and  frequency  m(k)  .  If  the 
fluid  is  compressible  so  that  it  sustains  both  pressure  (p)  and  density  (p) 

~  r,  ~ 

perturbations,  Dp/lt  =  c^  3p/3t,  the  perturbation  can  propagate,  e.g., 
in  form  of  longitudinal  sound  waves  with  finite  speed  =  (yPq/p  ) 2  and 
dispersion  law 

„ ?  =  cz  k*  +  i(5>'/3p  )mk? 
s  o 

In  a  fluid  with  a  viscosity  a  perturbation  may  also  propagate  in  form 
of  a  t i ansverse  or  shear  wave.  If  one  applies  the  curl  operation  to  the 
im  1'inpi  ess  i  hie  or  compressible  Naviet  -F>l"kcs  equation,  a  dispersion  law  is 
to  Mid  tor  the  transverse  perturbations  whieb  does  not  repieseut  a  wave 
phenomenon  hut  an  aperiodie  damping  process  witii 


=  -  (./p  )k- 
o 


As  is  known,  the  acoustic  dispersion  law  is  derived  from  a  hyperbolic  wave 
eipiation,  whereas  the  damping  constant  lot  t ho  transverse  modes  follows 
from  the  parabolic  vorticity  equation  (which  is  the  same  for  incompressible 
and  compressible  fluids).  From  experiments,  however,  it  is  established  that 
transverse  perturbations  (V  x  v^  =  ik  x  v^  i  tf)  propagate  as  (hyperbolic) 


HQ 


shear  waves  with  finite  speed. 

It  will  be  demonstrated  herein  that  the  (incompressible  or  compressible) 
Navier-Stok.es  equations  do  not  provide  a  correct  description  of  shear  waves. 

The  discrepancy  between  the  Navier-Stokes  equations  and  the  experiments  on 
shear  waves  is  resolved  by  introducing  viscous  stress  relaxation,  which  leads 
to  a  hyperbolic  transport  equation  for  shear  waves  in  incompressible  or  com¬ 
pressible  fluids.  For  this  reason,  the  transverse  or  shear  waves  represent 
"stress  relaxation  waves". 

As  an  illustration,  two  hyperbolic  ini tial -boundary-value  problems  for 

shear  waves  with  stress  relaxation  are  solved.  The  solutions  of  the  generalized 

Navier-Stokes  equations  with  stress  relaxation  represent  transverse  waves 

which  are  discontinuous  at  the  wave  front  and  have  a  finite  wave  speed, 

c  =  (p/p  t)'4  <  <«  (t  is  the  stress  relaxation  time).  The  first  treats  the 
o 

propagation  of  a  shear  wave  into  a  semi-infinite  fluid  space,  x  -  0,  produced 
by  a  temporal  velocity  impulse  at  the  boundary  x  -  0  (accelerated  wall).  The 
second  is  concerned  with  the  propagation  of  n  shear  wave  into  an  infinite 
fluid  space,  -  m  -  x  -  +  caused  by  a  spatial  velocity  pulse  in  the  plane 
x  =  0  at  time  t  =  0.  Both  solutions  are  valid  for  nonlinear  shear  waves  of 
arbitrary  intensity. 


i 
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PHYSICAL  PRINCIPLES 

In  conventional  fluid  mechanics,  — ^  it  is  assumed  that  inhomogeneities 
V  v  in  the  velocity  components  v.  produce  instantaneously  viscous  stresses 

i-  .1  i 

1!  .  Mathematically,  this  is  expressed  through  the  phenomenological  "flux"  ~ 

i  j 

"force"  relation 

II..  =  -  ,(V .  v .  +  V.v.  -  ~  V  v  6.  .) 
ij  i  j  j  x  3  k  k  i.| 

f- 

where  !•  is  the  viscosity  and  6  is  the  unit  tensor.  In  a  real  continuum, 
however,  velocity  inhomogeneities  do  not  switch  on  viscous  stresses  instan¬ 
taneously  but  in  accordance  with  a  relaxation  process  of  characteristic  time 
t.  By  means  of  the  kinetic  theory  of  gases  -\nnd  liquids,  --^one  can  show 
that  the  transport  equation  for  the  viscous  stresses  has  the  form  of  a 
temporal  (3/Dt)  and  spatial  (v-V>  relaxation  equation. 


-  -  II .  .  +  v,  V,  II  .  . 
Jt  J.J  k  k  jj 


t  Si..-  I  i  '(V.v.  +  V.v.  -  '  V  v  'S  .  . )  . 
ij  l  j  j  j  i  k  k  ij 


This  equation  is  approximate  insofar  as  the  the  coupiing  oi  heat  flows  <|^ 

and  stresses  II..  and  higher  order  terms  in  the  derivatives  of  v.  are 
ij  i 

2  3) 

neglected.— It  has  temporal  and  spatial  derivatives  as  required  for  a 
J-L  dependent  field  equation  and  is  Cali  lei  covariant.  If  relaxation  effects 
are  disregarded,  it  reduces  to  the  static  stress  equation. 

Thus,  consideration  ol  viscous  stress  relaxation  lends  to  a  refornmlat I  >n 
of  the  conventional  Nav fer-Stokes  theory  of  incompressible  and  compressible 
I luids.  In  place  of  the  Navi er-Stokes  equations,  we  have  the  hydrodynamic 
equations  with  viscous  stress  relaxation: 

p(ft  +  ''•VV')  ~  V>'  '  •  (1) 


—  +  v.Vfi  =  -  pV-v 

lit 


an> .  ->  4-y ,  ‘ii 


+  v-vn  + 


-  =  -  -  (Vv  +  Vv  - 
x  x  ' 


Eqs.  (1)  -  (3)  hold  for  incompressible  (V*v  -  0)  and  compressible  (  V  •  v  j-  0) 

fit  .  For  noniso thermal  systems,  the  transport  equations  for  thermal  energy 

A  51 

and  heat  flux  have  to  be  added  to  Eqs.  (1)  -  (3).  — — 

If  p  and t  can  be  treated  as  r-independent ,  It  is  mathematically  more 
convenient  to  use  instead  of  the  tensor  equation  (3)  the  vector  equation, 

—  v-n  +  v*  (v-vn)  +  t  v-n  =  ~y.i  (v2v  +  ~  w*v) ,  (A) 

since  Eq .  (1)  contains  the  force  density  V-TT.  If  temporal  and  spatial  re¬ 
laxation  of  the  viscous  stresses  is  disregarded,  Eqs.  (1)  and  (A)  combine  to 
the  classical  Navier-S tokes  equation,  p(3v/St  +  v*Vv)  =  -  7p  +  pV2v  +  ( lj./3)VV* 
Equations  (1)  -  (3)  represent  a  hyperbolic  system  both  in  the  compressible 
and  incompressible  cases.  On  the  other  hand,  the  conventional  incompressible 
Navier-Stokes  equations  are  parabolic.  The  corresponding  field  equations  for 
incompressible  fluids  are  obtained  by  setting  V-v  :  0  in  Eqs.  (2),  (3),  and  (A) 


t  > 


9:> 


INITIAL-BOUNDARY-VALUE  PROBLEM  FOR  v (0 , I ) -PULSE 


A  simple  method  for  the  generation  of  transverse  waves  in  a  viscous  fluid 

consists  in  setting  the  plane  x  =  0  bounding  a  semi-infinite  fluid  (x  _  0, 

|  s'  1  v  >>■,  I  z  I  1  "•)  into  sudden  motion  v  =  v  H(t)e  ,  where  H(t)  is  the 
1,11  w  y 

Heaviside  step  function.  The  resulting  viscous  interaction  between  the  fluid 

and  the  accelerated  wall  produces  a  curl  n  x  [y]  -  e  v(x  =  Q,t)  at  the  fluid 

surface  which  propagates  in  form  of  a  transverse  wave  through  the  fluid  in 

the  x-direction.  In  this  dynamic  process,  the  fluid  velocity  is  of  the  form 

v  =  10,  v(x,t),  0}  so  that  V • v  =  3v/3 y  =  0  and  v*Vv  =  0,  i.e.,  the  fluid 

motion  behaves  incompressible  (even  if  the  fluid  is  compressible)  and  linear. 

->  ■<  — ►  — >■  ^ — > 

Furthermore,  vVJI  =  v  DR/3y  =  0  since  )l  has  only  a  single  component  =  H(x,t) 
by  Eq.  (3),  and  Vp  =  (5  by  Eq.  (1). 

Thus,  Eqs.  (1)  -  (3)  lead  to  the  following  initial-boundary-value  problem 

-V 

for  the  transverse  velocity  wave  v(x,t)  in  the  y-direction  propagating  in 
the  x-direction,  as  a  result  of  the  sudden  wall  motion  in  the  plane  x  =  0: 


3  v  _311_ 

Po  It  dx 

(5) 

on  n  _  _  p 

0 1  t  'i  3  x 

♦ 

<&) 

v(x  =  U,t)  =  v  H(t),  t  1  0 

> 

(7) 

o 

A 

o 

1! 

/**N 

o 

II 

> 

> 

(8) 

>  v  ( x ,  t  =  0)/3t  =  0,  x  •  0 

(9) 

where  H(t)  =0,  L  -  -  0,  and  H(t)  =  i,  t  -  +  0. 

Equations  (5)  -  (6)  represent  a  hyperbolic  system,  from  which  one  obtains 
by  elimination  wave  equations  for  the  stress  component  11  E  n(x,t)  and  the 

*y 

velocity  lield  v(x,t): 


3  7 II 


=  c?  - 
fix 


;v  ri 

ot7 


1  Oli 

t  3t 


(10) 


whore 


C  =  (;,/p  T)2  (12) 

o 

is  the  (maximum)  spued  of  the  stress  relaxation  wave.  Both  i.  (x,t)  and  v(x,t) 
satisfy  similar  (hyperbolic)  wave  equations  with  the  same  wave  speed  c.  In 
the  limit,  t  ->  0  and  c  -►  ",  with  c?t  ►  i  /(  ,  Eqs.  (10)  and  (11)  reduce  to 

parabolic  equations,  according  to  which  boundary  values  of  Il(x,t)  and  v(x,t) 
would  diffuse  with  infinite  speed  into  the  fluid  (conventional  Navi er-Stok.es 


theory).  Accordingly,  only  for  t  >  0  and 


■,  transverse  or  shear  waves 


exist  in  the  fluid  which  represent,  therefore,  stress  relaxation  waves. 

According  to  Eq .  (11)  and  Eqs.  (7)  -  (9),  the  velocity  field 
v(x,t)  ~  v  u(f, t)  of  the  stress  relaxation  wave  under  consideration  is  des¬ 
cribed  by  the  dimensionless  initial-boundary-vnlue  problem: 

32u  3u  _  32u 

+  3*  ”  al7  ,  (13) 

u(4  =  0,  t)  =  H(t),  t  0,  ,  (14) 

u ( f ,  t  =  0)  =  0,  f,  >  0  ,  (15) 

3u(C,  £  =  0)/D t  =  0,  4  >  0  ,  (16) 


where 

U (\,t)  =  v(x,t)/v,  4  =  x/cr,  t  =  t/t 

Equation  (13)  -  (16)  are  solved  by  means  of  the  Laplace  transform 

i  •  6) 

technique  —  winch  gives 

u(4,s)  -■  L(u(f.,f)]  =  /  e  s!"  u(4,()d( 

0 

u((),s)  =  L [u(n,^)  ]  =  J  e  H(Od(  =  s 


F/6  20/9 
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d2u 

w 


Since  the  initial  conditions  (15)  -  (16)  vanish,  Eqs.  (13)  -  (14)  yield  for 
the  transformed  velocity  u(£,s)  the  ordinary  boundary-value  problem, 

(s2  +  s)u  =  0  ,  (20) 

u(£  =  0,s)  =  s  ^  .  (21) 

Since  u(£,s)  must  be  finite  for  £  -*■  “,  the  solution  of  Eqs.  (20)  -  (21)  is 

(22) 

The  inverse  Laplace  transform  gives  for  the  velocity  field  the  complex 
integral  , . 

°  7+1“  , 

_ 1  _  / e2  i  2  .  e/* 

+  s;  f’eSlds.  (23) 


-  .  .  -1  -(s2  +  s)^  r 

u(£,s)  =  s  e 


u (£,t)  =  (1/2*1)  J  s  e  ^ 

Y~i“ 

Hence, 

u  (£,*)  =  3$(£,-t)/3£ 

where 

y+ioo 

<Kf,,'t)  =  -  2~r  /  s'1  f  (£,s)eS*ds 


Y-i00 


(24) 


(25) 


and 


f(£,s)  =  e  (S2+S)W  +  a)*  .  (26) 

According  to  a  known  inversion  integral,—^  the  inverse  transform  of  Eq.  (26)  is 

(27) 


fU,L)  =  r1lf(C,s)]  =  lo(h (t2  -  £2)S  H(t  -  £) 


where  I  U)  is  the  modified  Bessel  function  of  order  v.  By  Eq  (25), 


<!>(£,£)  =  -  L  1  [s  1  f(£,s)]  =  -  J  L  X[f(£,s)]di 

t  ° 

=  -  /  f  (£  ,-£)dC 
0 


i.e. 


♦(£,*)  =  -  H u  -  o  J  e'^^I  (^(a2  -  C2)  *)da 


(28) 


(29) 


From  this  potential,  the  dimensionless  velocity  field  is  obtained  as 

u(£ ,t)  =  H(*  -  C)  [e”^  *  +  hK  J  e~Hu  (a2  -  Z2)~H  I,(>s(o2  -  52)*‘)do]  (30) 

£ 
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in  accordance  with  Eq.  (24).  By  Eq.  (17),  the  corresponding  dimensional 
solution  for  the  velocity  field  is 
v(x,u)  =  vH(cl  -  x)  {e 

t/ 1  ,  ,  i 

+  (x/2ct)  /  e"'?U  [cc2-  (x/ct)2]  2  I^lu2  -  (x/ct)2]  2)da j.  (31) 

x/cx 

Eq.  (31)  indicates  that  the  transverse  stress  relaxation  wave  is  dis¬ 
continuous  at  x  =  ct,  *-he  position  of  the  wave  front  at  time  t.  At  any  time 
0  i  t  S,  only  the  region  0  1  x  1  ct  of  the  fluid  is  excited  by  the  wave, 
since  v(x  <  ct,  t)  >0  and  v(x  >  ct,  t)  =  0  by  Eq.  (31).  The  velocity  signal 
v(0,t)  =  v  H(t)  generated  at  the  boundary  x  =  0  at  time  t  is  thus  transported 
with  finite  speed  c  =  (p  /pot) 2  <  °°  in  form  of  a  discontinuous  wave  into  the 
fluid  space  x  >  0  as  t  >  0  increases. 


2 

Application  of  the  asymptotic  formula  1  (z)  ~  e  /(2:rz)2,  |z|  >>  1,  and 
expansion  of  2  =  ^[a2  -  (x/ct)2]2  for  large  ;x-values  in  Eq.  (3.1)  yields,  in 
the  limit  t  ->  0,  ct  -*  0: 

00  2  1 

v(x,t)  =  v(2//tT)  /  e  6  dE,  n  =  x/2(—  t)'2  .  (32) 

n  ‘o 

This  is  the  familiar  solution  of  the  parabolic  Navier-Stok.es  equations  due 

to  Stokes.  Eq.  (32)  suggests  that  v(x,t)  >  0  throughout  the  entire  fluid 

0  •  x  *  for  any,  no  matter  how  small  time  t  >  0.  Thus,  the  parabolic 

Stokes  solution  gives  a  completely  misleading  picture  for  a  shear  wave  in 

form  of  a  diffusion  process  which  spreads  with  infinite  speed. 

0  ]  2 

Fig.  1  shows  u(c,t)  versus  £  for  t  =  10  ,  10  ,  10  ,  with  wave  fronts  at 
t  =  10  ,  10  ,  10~.  It  is  seen  how  the  perturbation  u(0,i)  =  H(f)  produced  at  the 
wall  -  0  moves  in  form  of  a  discontinuous  wave  into  the  fluid  space  Z  -  0 
so  that  an  increasing  but  finite  region  0  1  (,  -  r,  of  the  fluid  is  set  into  motion 

with  increasing  t.  in  tile  limit  (  =  u (f.,0  =  1  throughout  the  fluid, 

<  < 

j 


0 
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INITIAL-BOUNDARY-VALUE  PROBLEM  FOR  v(x,0) -PULSE 


Another  fundamental  method  for  shear  wave  generation  makes  use  of  a  velocity 


pulse  vo=v(x,0)e^  generated  of  time  t  =  0  within  a  limited  region  |x|  <  Ax. 
The  decay  of  this  velocity  pulse  occurs  in  form  of  a  shear  wave  with  velocity 
field  v  =  {0,  v(x,t),  0}  in  the  y-direction  propagating  in  the  x-directions . 

f  -►  -V  •>  •< — f  4  <■-> 

Accordingly,  V •  v  =  3v/3y  =  0,  v-Vv  =  0,  and  vVIl  =v3Ii/Dy  =  0,  since  II  has 

only  a  single  component  II  =  n(x,t).  Again,  the  stransverse  wove  "behaves" 

xy 

incompressible  and  linear,  and  Vp  =  3  by  Eq.  (1). 

As  in  the  previous  problem,  Eqs.  (1)  -  (3)  give  the  wave  Eqs.  (10)  and 
(11)  for  n(x,t)  and  v(x,t),  respectively.  Hence,  the  shear  wave  produced  by 
the  velocity  pulse  v  =  v(x,0)e  is  described  by  the  initial-boundary-value 

a  y 


problem: 


3  2v  1  3  v  2 
Ft2  T  Ft  =  c  a? 


v ( x , t  =  0)  =  VQ(x) 


I  v ( x ,  t  =  0)  /U  t  =  wQ(*)  ,  M 


< 

—  CO 


(33) 

(35) 

(35) 


where  |w  (x) I  -  0  is  included  for  reasons  of  generality.  The  solution  of 
o 

g) 

Eqs.  (33)  -  (35)  is  accomplished  by  means  of  Riemann's  method,  — 

_  x+ct 
v  ( x ,  t )  =  e  T  (i;[vo(x-ct)  +  v^(x+ct)]  +  %  /  Y(x,  t,a)da}  (36) 

x-ct 

where 


T(x,  t ,  ct )  =  vQ(a)  (t/2t)  lj_(2“  [c2t2  -  (u  -  x)2  ]  2)  /  [c2t2  -  (a  -  x)2]2 
+  -  [w  (a)  +  tv-1-  v  (a)]  I  (^-  [c2t2  -  (a  -  x)2]2)  .  (37) 

CO  o  o  ZCT 

As  a  concrete  example  tor  the  initial  conditions  in  Eqs.  (34)  -  (35),  an 
initial  velocity  distribution  of  the  form  of  a  Dirac  pulse  is  chosen, 

v  (x)  =  v  6(x),  w  (x)  =  0,  ™  .  (38) 

o  o  o 
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In  this  case,  the  general  solution  in  Eqs.  (36)  -  (37)  becomes,  in  dimen¬ 
sionless  form, 

u U,t)  =  %  e~l  [6(5-  t)+  6(5+  t) 

+  Ijilt?-  e,2]h+lQ([t2  -  52]4)],  1 5 1  1 

=  0  ,  Id  '  *,  (39) 

where 

u(5,t)  =  v(x,t)/(vo/2ci) ,  5  =  x/2cx,  t=t/2i  .  (40) 

Eq.  (39)  indicates  that  the  shear  wave  spreads  in  the  space  | 5 |  i  ■»  in 
form  of  a  symmetrical  wave,  u(-  5 ,-£)  =  u(+  E,,t)  due  to  the  symmetry  of  the 
initial  conditions  (38).  The  wave  is  discontinuous  at  its  fronts  5  =  ±  t, 
which  propagate  with  the  speed 

v (x , t )  =  ±  c,  x  =  ±  ct  .  (41) 

In  the  limit  c  ■+  0,  ct  -►  0,  application  of  the  asymptotic  formula 

2  J,  L 

I^(z)  ~  e  /  (2nz)  2,  |z|  >>  1,  and  expansion  of  z  =  [a2  -  52]2  for  large  '(-values 
in  Eq.  (39)  yields 

u(5  »-0  =  (2irf)-%  1 5 1  -  -  .  (42) 

This  is  the  corresponding  solution  of  the  parabolic  Navier-Stokes  equations. 

Eq.  (42)  would  indicate  that  the  shear  wave  has  the  form  of  a  Gaussian  ex¬ 
tending  from  5  =  _0°  to  5=  +  °°  for  any,  no  matter  how  small  time  t  >  0  (cor¬ 
responding  to  an  infinite  speed  of  propagation).  It  is  obvious  that  the 
solution  (42)  is  physically  not  meaningful. 

In  Fig.  2,  the  dimensionless  velocity  field  u(5,-t)  of  the  shear  wave  is 
shown  versus  5  for  t  =  10®,  10^,  102 ,  the  wave  fronts  being  in  each  case  at 
5=1  t.  Due  to  the  finite  wave  speed  c,  the  fluid  is  not  excited  in  the 
region  |f|  >  f  ahead  of  the  wave  fronts.  The  shape  of  the  wave  is  flat 


% 
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with  relatively  steep  flanks  leading  to  the  discontinuous  fronts.  Thus, 
the  shear  wave  does  not  resemble  the  Gaussian  of  the  parabolic  theory,  Eq. 


i 

1 


i 


(40). 
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CONCLUSIONS 

A  generalization  of  the  Navier-Stokes  equations  is  presented  considering 
viscous  stress  relaxation,  which  results  in  a  physically  meaningful  theory  for 
transverse  waves  in  viscous  fluids.  The  fundamental  speed  of  the  stress 
relaxation  waves  is  given  by  c  =  (p/px)J,  where  \i  is  the  viscosity,  p  is  the 

9) 

density,  and  t  is  the  relaxation  time  of  the  stress  tensor.  For  any  medium, 

it  is  c  ~c  ,  where  c  is  the  speed  of  the  longitudinal  waves,  e.g., 

c  =  1.3  x  10 cm  sec  ^  and  c  *  1.5  x  10^  cm  sec  ^  for  water  at  T  =  20°C  and 

s 

P  =  1  atm. 
o 

Exact  solutions  are  derived  for  stress  relaxation  waves  propagating  in 
the  x-direction  due  to  velocity  pulses  v(0,t)  and  v(x,0)  in  the  y-directlon, 
respectively.  For  the  geometry  of  these  transverse  waves,  the  nonlinear, 
generalized  Navier-Stokes  equations  become  linear,  so  that  the  solutions  given 
hold  for  waves  of  arbitrary  intensity.  The  solutions  are  discontinuous  at 
the  wave  fronts,  which  is  typical  for  hyperbolic  field  equations.  The  corres¬ 
ponding  solutions  of  the  conventional  Navier-Stokes  equations  indicate  a 
diffusion  process  with  infinite  wave  speed,  i.e.,  give  a  qualitatively  and 
quantitatively  insufficient  picture  of  the  propagation  of  transverse  waves 
in  fluids. 

In  the  simplified  stress  relaxation  equation  (3)  proposed,  the  term 

'[•Vv  neglected  since  it  is  of  the  order-of-magnitude  of  (h/t)jVv|2, 

which  is  nonlinear  in  the  derivatives.  It  should  be  noted  that  the  term 

H*Vv  vanishes  exactly  for  the  wave  problems  treated  above,  Il*Vv=0  since  v  = 

{0,v(x,t),0}  and  II  has  only  a  single  component  n  (x,t).  For  this 

xy 

reason,  the  solutions  presented  arc  exact  solutions  of  the  nonlinear 
Navier-Stokes  equations  with  viscous  stress  relaxation. 
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